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Abstract 

We investigate a coupled system of the Ricci flow on a closed manifold M with the 
harmonic map flow of a map cj> from M to some closed target manifold N, 

§- t g = -2Rc + 2aV0 ® V(j>, §- t 4> = T g <j>, 

where a is a (possibly time-dependent) positive coupling constant. Surprisingly, the 
coupled system may be less singular than the Ricci flow or the harmonic map flow alone. 
In particular, we can always rule out energy concentration of cj> a-priori by choosing a 
large enough. Moreover, it suffices to bound the curvature of (M,g(t)) to also obtain 
control of <f> and all its derivatives if a > a > 0. Besides these new phenomena, the 
flow shares many good properties with the Ricci flow. In particular, we can derive the 
monotonicity of an energy, an entropy and a reduced volume functional. We then apply 
these monotonicity results to rule out non-trivial breathers and geometric collapsing at 
finite times. 

Resume 

Nous etudions un systeme d'equations consistant en un couplage entre le flot de Ricci 
et le flot harmonique d'une fonction <f> allant de M dans une variete cible N, 

§- t g = -2Rc + 2qV0 ® V(j>, §- t< j> = T g <f>, 

ou a est une constante de couplage strictement positive (et pouvant dependre du 
temps). De maniere surprenante, ce systeme couple peut etre moins singulier que 
le flot de Ricci ou le flot harmonique si ceux-ci sont consideres de maniere isolee. En 
particulier, on peut toujours montrer que la fonction cj> ne se concentre pas le long de 
ce systeme a condition de prendre a assez grand. De plus, il est suffisant de borner 
la courbure de (M,g(t)) le long du flot pour obtenir le controle de cj> et de toutes ses 
derivees si a > a > 0. A part ces phenomenes nouveaux, ce flot possede certaines 
proprietes analogues a celles du flot de Ricci. En particulier, il est possible de montrer 
la monotonie d'une energie, d'une entropie et d'une fonctionnelle volume reduit. On 
utilise la monotonie de ces quantites pour montrer l'absence de solutions en "accordeon" 
et l'absence d'effondrement en temps fini le long du flot. 
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1 Introduction and main results 



Let (M m ,g) and (N n , 7) be smooth Riemannian manifolds without boundary. According 
to Nash's embedding theorem [30] we can assume that N is isometrically embedded into 
Euclidean space (N n ,j) °-t> M. d for a sufficiently large d. If ■ N — s- R d denotes this 
embedding, we identify maps : M N with ejv '• M — > R d , such maps may thus be 
written as — (4> x )i<\<d- Harmonic maps : M N are critical points of the energy 
functional 

E(0) = f |V0| 2 dV. (1.1) 

JM 

Here, |V0| 2 := 2e(0) = g ij V 'j0 A V j-0 A denotes the local energy density, where we use the 
convention that repeated Latin indices are summed over from 1 to to and repeated Greek 
indices are summed over from 1 to d. We often drop the summation indices for when 
clear from the context. Harmonic maps generalize the concept of harmonic functions and in 
particular include closed geodesies and minimal surfaces. 

To study the existence of a harmonic map homotopic to a given map 0o : M — > N, Eclls 
and Sampson [13] proposed to study the L 2 -gradient flow of the energy functional (1.1), 

£0 = T ff 0, 0(0) = 0o, (1-2) 

where r ff denotes the intrinsic Laplacian of 0, often called the tension field of 0. They 
proved that if N has non-positive sectional curvature there always exists a unique, global, 
smooth solution of (1.2) which converges smoothly to a harmonic map 0oo : M — > N homo- 
topic to 0o as t — y 00 suitably. On the other hand, without an assumption on the curvature 
of N, the solution might blow up in finite or infinite time. Comprehensive surveys about 
harmonic maps and the harmonic map flow are given in Eclls-Lcmaire [11, 12], Jost [18] and 
Struwe [40]. The harmonic map flow was the first appearance of a nonlinear heat flow in 
Riemannian geometry. Today, geometric heat flows have become an intensely studied topic 
in geometric analysis. 

Another fundamental problem in differential geometry is to find canonical metrics on Rie- 
mannian manifolds, for example metrics with constant curvature in some sense. Using the 
idea of evolving an object to such an ideal state by a nonlinear heat flow, Richard Hamilton 
[14] introduced the Ricci flow in 1982. His idea was to smooth out irregularities of the 
curvature by evolving a given Riemannian metric g on a manifold M with respect to the 
nonlinear weakly parabolic equation 

f.g = -2Rc, g(Q) = g , (1-3) 

where Rc denotes the Ricci curvature of (M,g). Strictly speaking, the Ricci flow is not the 
gradient flow of a functional J(g) = J M F(d 2 g, dg, g)dV, but in 2002, Perelman [31] showed 
that it is gradient-like nevertheless. He presented a new functional which may be regarded 
as an improved version of the Einstein-Hilbert functional E(g) = J M RdV, namely 

3(9, f)~ J^(R+\Vf\ 2 )e-fdV. (1.4) 

The Ricci flow can be interpreted as the gradient flow of ? modulo a pull-back by a fam- 
ily of diffeomorphisms. Hamilton's Ricci flow has a successful history. Most importanty, 
Perelman's work [31, 32] led to a completion of Hamilton's program [15] and a complete 
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proof of Thurston's geometrization conjecture [42] and (using a finite extinction result from 
Perelman [33] or Colding and Minicozzi [8, 9]) of the Poincare conjecture [34]. Introductory 
surveys on the Ricci flow and Pcrelman's functionals can be found in the books by Chow and 
Knopf [4], Chow, Lu and Ni [7], Miiller [27] and Topping [43]. More advanced explanations 
of Perelman's proof of the two conjectures are given in Cao and Zhu [3] Chow et al. [5, 6], 
Kleiner and Lott [19] and Morgan and Tian [25, 26]. A good survey on Perelman's work is 
also given in Tao [41]. 

The goal of this article is to study a coupled system of the two flows (1-2) and (1.3). Again, 
we let (M m ,g) and (TV™, 7) be smooth manifolds without boundary and with (TV™, 7) °->- M. d . 
Throughout this article, we will assume in addition that M and TV are compact, hence 
closed. However, many of our results hold for more general manifolds. 

Let g(t) be a family of Riemannian metrics on M and <p(i) a family of smooth maps from M 
to TV. We call (g(t), 4>(t))te[o,T) a solution to the coupled system of Ricci flow and harmonic 
map heat flow with coupling constant a(t), the (RH) a flow for short, if it satisfies 



Here, r g <p denotes the tension field of the map <p with respect to the evolving metric g, and 
a(t) > denotes a (time-dependent) coupling constant. Finally, V0 (8 V</> has the com- 
ponents (V(j> <8> V0)tj = V;(T> A Vj</> A . In particular, |V0| 2 as defined above is the trace of 
® with respect to g. 

The special case where TV C R and a = 2 was studied by List [22], his motivation coming 
from general relativity and the study of Einstein vacuum equations. Moreover, List's flow 
also arises as the Ricci flow of a warped product, see [28, Lemma A. 3]. After completion of 
this work, we learned that another special case of (RH) a with TV C SL(kM)/SO(k) arises 
in the study of the long-time behaviour of certain Type III Ricci flows, see Lott [23] and a 
recent paper of Williams [44] for details and explicit examples. 

The paper is organized as follows. In order to get a feeling for the flow, we first study 
explicit examples of solutions of {RH) a as well as soliton solutions which are generalized 
fixed points modulo diffcomorphisms and scaling. The stationary solutions of (RH) a satisfy 
Rc = aV0 <S> V</>, where 4> is a harmonic map. To prevent (M,g(t)) from shrinking to a 
point or blowing up, it is convenient to introduce a volume-preserving version of the flow. 

In Section 3, we prove that for constant coupling functions a(t) = a > the (RH) a flow can 
be interpreted as a gradient flow for an energy functional I?a(ff, 4>, /) modified by a family of 
diffcomorphisms generated by V/. If (g(t),<p(t)) solves (RH) a and e~f is a solution to the 
adjoint heat equation under the flow, then 7 a is non-decreasing and constant if and only 
if (g(i), 4>(t)) is a steady gradient soliton. In the more general case where a(t) is a positive 
function, the monotonicity result still holds whenever a(t) is non- increasing. This section is 
based on techniques of Perelman [31, Section 1] for the Ricci flow. 

In the fourth section, we prove short-time existence for the flow using again a method from 
Ricci flow theory known as DeTurck's trick (cf. [10]), i.e. we transform the weakly parabolic 
system (RH) a into a strictly parabolic one by pushing it forward with a family of diffcomor- 
phisms. Moreover, we compute the evolution equations for the Ricci and scalar curvature, 
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the gradient of <j) and combinations thereof. In particular, the evolution equations for the 
symmetric tensor Sij := Rij — and its trace S = R — a|V</>| 2 will be very useful. 

In Section 5, we study first consequences of the evolution equations for the existence or 
non-existence of certain types of singularities. Using the maximum principle, we show that 
min^gM S(x, t) is non-decreasing along the flow. This has the rather surprising consequence 
that if \X7cf)\ 2 (xk,tk) — > oo for % /• T, then R(xk,tk) blows up as well, i.e. g{tk) must 
become singular as tk /* T. Conversely, if |Rm| stays bounded along the flow, |V0| 2 must 
stay bounded, too. This leads to the conjecture that a uniform Ricmann-bound is enough 
to conclude long-time existence. This conjecture is proved in Section 6. To this end, we 
first compute estimates for the Riemannian curvature tensor, its derivatives and the higher 
derivatives of <j> and then follow Bando's [1] and Shi's [38] results for the Ricci flow to derive 
interior-in-timc gradient estimates. 

In Section 7, we introduce an entropy functional W Q (<7, 0, /, r) which corresponds to Percl- 
man's shrinker entropy for the Ricci flow [31, Section 3]. Here r = T — t denotes a backwards 
time. For a(t) = a > 0, the entropy functional is non-decreasing and constant exactly on 
shrinking solitons. Again, the entropy is monotone if we allow non-increasing positive cou- 
pling functions a(t) instead of constant ones. Using 3 a and W Q we can exclude nontrivial 
breathers, i.e. we show that a breather has to be a gradient soliton. In the case of a steady 
or expanding breather the result is even stronger, namely we can show that <f>(t) has to be 
harmonic in these cases for all t. 

Finally in the last section, we state the monotonicity of a backwards reduced volume quantity 
for the (RH) a flow with positive non- increasing a(t). This follows from our more general 
result from [29]. We apply this monotonicity to deduce a local non-collapsing theorem. 

In the appendix, we collect the commutator identities on bundles like T*M <E> <p*TN, which 
we need for the evolution equations in Section 4 and 6. 

This article originates from the authors PhD thesis [28] from 2009, where some of the proofs 
and computations are carried out in more details. The author likes to thank Klaus Ecker, 
Robert Haslhofer, Gerhard Huiskcn, Tom Ilmanen, Peter Topping and in particular Michael 
Struwe for stimulating discussions and valuable remarks and suggestions while studying 
this new flow. Moreover, he thanks the Swiss National Science Foundation that partially 
supported his research and Zindine Djadli who translated the abstract into flawless French. 

2 Examples and special solutions 

In this section, we only consider time-independent coupling constants a(t) = a. First, we 
study two very simple homogeneous examples for the (RH) a flow system to illustrate the 
different behavior of the flow for different coupling constants a. In particular, the existence 
or non-existence of singularities will depend on the choice of a. We study the volume- 
preserving version of the flow as well. We say that (g(t), 4>{t)) is a solution of the normalized 
(RH) a flow, if it satisfies 

[ ^. g = -2Rc + 2«V0®V ( /»-f^ 5 / (R-a^^dV, 

< Jm (2.1) 

l jl<t> = Tg<t>. 
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Two homogeneous examples with cf) = id 

Assume that (M,g(0)) is a round two-sphere of constant Gauss curvature 1. Under the 
Ricci flow, the sphere shrinks to a point in finite time. Let us now consider the (RH) a 
flow, assuming that (N, 7) = (M,g(Q)) and </>(0) is the identity map. With the ansatz 
g(t) = c(t)g(0), c(0) = f and the fact that (f>(t) = (f)(0) is harmonic for all g(t), the (RH) a 
flow reduces to 

■§- t c{t) = -2 + 2a. 

For a < 1, c(t) goes to zero in finite time, i.e. (M, g(t)) shrinks to a point, while the scalar 
curvature R and the energy density |V</>| 2 both go to infinity. For a = I, the solution is 
stationary. For a > 1 , c(t) grows linearly and the flow exists forever, while both the scalar 
curvature R and the energy density |V0| 2 vanish asymptotically. Instead of changing a, we 
can also scale the metric 7 on the target manifold. Mapping into a larger sphere has the 
same consequences as choosing a larger a. Note that the volume-preserving version (2.1) of 
the flow is always stationary, as it is for the normalized Ricci flow, too. 

A more interesting example is obtained if we let (M 4 ,g(t)) = (S 2 x L,c(t)g§2 © d(t)g^), 
where (§ 2 ,g§2) is again a round sphere with Gauss curvature 1 and (L,g^) is a surface (of 
genus > 2) with constant Gauss curvature —1. Under the Ricci flow, ^c(t) = —2 and 
■§[d(t) = +2. In particular, c(t) goes to zero in finite time while d(t) always expands. Under 
the normalized Ricci flow J^g = — 2Rc + ^ g j- R dV, we have 

£c = -2 + ^ = -f-c 2 , JL rf = + 2 + ^ = + I + rf 2 . 

Again, c(t) goes to zero in finite time. At the same time, d(t) goes to infinity. Now, let us 
consider the (RH) a flow for this example, setting (AT, 7) = (M, g(0)) and (f)(0) = id. First, 
note that 0(0) is always harmonic and thus <f>(t) = 0(0) is unchanged. The identity map 
between the same manifold with two different metrics is not necessarily harmonic in general, 
but here it is. The flow equations reduce to J^c(i) = — 2 + 2a and §id(t) = +2 + 2a. While 
d(t) always grows, the behavior of c(t) is exactly the same as in the first example above, 
where we only had a two-sphere. On the other hand, if we consider the normalized flow 
(2.1), we obtain 

= (a - 1) - (a + f)c 2 , = (a + 1) - (a - l)d 2 . 

In the case where a < 1, c(i) goes to zero in finite time, while d(t) blows up at the same 
time, similar to the normalized Ricci flow above. For a = I, c(t) = (1 + 2t)~ 1 goes to zero 
in infinite time while d(t) = (1 + 2t) grows linearly, i.e. we have long-time existence but no 
natural convergence (to a manifold with the same topology). In the third case, where a > I, 
both c(t) and d(t) converge with 

c(t) ^/g, d(t) -+ y/g, as t -J- 00. 

These examples show that both the unnormalizcd and the normalized version of our flow 
can behave very differently from the Ricci flow if a is chosen large. In particular, they may 
be more regular in special situations. 

Volume-preserving version of (RH) a 

Here, we show that the unnormalizcd (RH) a flow and the normalized version are related by 
rescaling the metric g and the time, while keeping the map (f> unchanged. Indeed, assume 
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that (g(t), 4>it)) t ^[o,T) is a solution of (RH) a . Define a family of rescaling factors X(t) by 

/ r s-2/m 

X{t):=lJ^dV g(t) \ , ie[0,T), (2.2) 

and let g(t) be the family of rescaled metrics g(t) — X(t)g(t) having constant unit volume 

/ dV m = I X m ' 2 {t) dV gW = 1, vt e [0,T). 
Jm Jai 

Write S = Rc - aV^> ® S7<j) with trace S = R- a|V0| 2 . It is an immediate consequence of 
the scaling behavior of Rc, R, V</> ® V</> and |V</>| 2 that § = § and S = A _1 5 for g = Xg. 
Note that X(t) is a smooth function of time, with 

2 + m 

j t x ^=-^{S M dv ^) m S A y s)dV9= ™ x fj dV9= ™ x2 fj dy v- 

Now, we rescale the time. Put t(t) := J* X(s)ds, so that ^| = X(t). Then, we obtain 

&9 = A- 1 ^ (Xg) = jk9 + (A~ 2 f X)g = -2$ + 2-gi S dV s , 

Jm 

J> = A -iJ> = A -i r ^ = T _^. 

This means that (g(t).<p(t)) solves the volume-preserving (RH) a flow (2.1) on [O.T), where 
T = Jo A(s)ds. 



Gradient solitons 

A solution to (RH) a which changes under a one-parameter family of diffcomorphisms on M 
and scaling is called a soliton (or a self-similar solution) . These solutions correspond to fixed 
points modulo diffcomorphisms and scaling. The more general class of periodic solutions 
modulo diffcomorphisms and scaling, the so-called breathers, will be defined (but also ruled 
out) in Section 7. 

Definition 2.1 

A solution (g(i), 0(i))te[o.T) °f {RH) a is called a soliton if there exists a one-parameter 
family of diffeomorphisms ip t '■ M —> M with ipo = idjvf and a scaling function c : [0, T) — > R+ 
such that 

f g{t) = c(t)V t *.9(0), 

The cases JjC = c < 0, c = and c > correspond to shrinking, steady and expanding 
solitons, respectively. If the diffeomorphisms ipt ore generated by a vector field X(t) that is 
the gradient of some function f(t) on M , then the soliton is called gradient soliton and f 
is called the potential of the soliton. 



Lemma 2.2 

Let {g(t), 4>{t))t£[a,T) be a gradient soliton with potential f. Then for any t$ £ [0,T"), this 
soliton satisfies the coupled elliptic system 

J = Rc - aV<j) <E)V<f) + Hess(/) +ag, 

1 o = t 5 0-(V0,v/>, (2 - 4) 



6 



for some constant o~(to). Conversely, given a function f on M and a solution of (2.4) at 
t = 0, there exist one-parameter families of constants c(t) and diffeomorphisms ip t : M — > AI 
such that defining (g(t),(f>(t)) as in (2.3) yields a solution of {RH) a . Moreover, c(t) can be 
chosen linear in t. 

Proof. Suppose we have a soliton solution to (RH) a . Without loss of generality, c(0) = 1 
and tpQ = idjvf • Thus, the solution satisfies 

-2Rc(. 9 (O))+2«(V0®V0)(O) = f<?(i)|t=o - §i (c(t)r t g(0))\t=o 

= c(0). 9 (0) + L x{0) g(0) = c(0)g(0) + 2 Hess(/(-, 0)), 

where X(t) is the family of vector fields generating ip t - Moreover, we compute 
M)(0) = I^)l*=o = ^x(o)0(O) = (V^VJ). 

Together, this proves (2.4) with a = |c(0) for to = 0. Hence, by a time-shifting argument, 
(2.4) must hold for any to £ [0,T). One can easily see that a(to) = c(to) / 2c(to) ■ 

Conversely, let (g(0), <^(0)) solve (2.4) for some function / on M. Define c(t) := 1 + 2at and 
X(t) := V//c(i). Let ipt be the diffeomorphisms generated by the family of vector fields 
X(t) (with Vo = idA/) and define (g(t),(j)(t)) as in (2.3). For a < this is possible on the 
time interval t £ [0, ^), in the case a > it is possible for t £ [0, oo). Then 

%g(t) = c(i)V?G?(0)) + c(t)^(L x{t) g(0)) = ^(2^(0) + £ v/5 (0)) 

= ip* t (2ag(Q) + 2Hcss(/)) = ipl (-2Rc( 5 (0)) + 2a(V0 ® V0)(O)) 
= -2 Rc(.g(i)) + 2a(V0 ® V0)(i), 

as well as 

&#t) = V t *(£x (t )0(O)) = #<V0(O), V//c(t)) - c(t)-Vt*(T fl( o)#0)) 

This means that (c/(i), </>(£)) is a solution of (RH) a and thus a soliton solution. □ 

By Lemma 2.2 and rescaling, we may assume that c(t) = T — t for shrinking solitons (here T 
is the maximal time of existence for the flow), c(t) = 1 for steady solitons and c(t) = t — T for 
expanding solitons (where T defines a birth time). An example for a soliton solution is the 
very first example from this section, where (M, g(0)) = (N, 7) = (S 2 , g§2 ) and </>(0) = id. For 
a < 1, the soliton is shrinking, for a = 1 steady and for a > 1 expanding. Since <pt = idjvf 
for all t in all three cases, these are gradient solitons with potential / = 0. 

Taking the trace of the first equation in (2.4), we see that a soliton must satisfy 

R — a|V0| 2 + A/ + am = 0. (2.5) 

Taking covariant derivatives in (2.4) and using the twice traced second Bianchi identity 
VjRij = h'ViR, we obtain (analogous to the corresponding equation for soliton solutions of 
the Ricci flow) 

R - a|V0| 2 + |V./f + 2a f = const. (2.6) 
Finally, with f(-,t) = V* (/(•, 0)), we get 

lf = Z x f=\Vf\ 2 . (2.7) 
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Combining this with (2.5), we obtain the evolution equation 

(f + A)/ = |V./f -R + q|V0| 2 - am. (2.8) 

For steady solitons, for which the formulas (2.4)-(2.8) hold with a = 0, equation (2.8) is 
equivalent to u = e~? solving the adjoint heat equation 

n*u = -§- t u - Au + Ru - a\W(f>\ 2 u = 0. (2.9) 

For shrinking solitons, (2.4)-(2.8) hold with a(t) = —\(T - i) -1 and (2.8) is equivalent to 
u = (4n(T — t))~ m / 2 e~f solving the adjoint heat equation. Finally, for expanding solitons, 
<?{t) = +l(t - T)- 1 and (2.8) is equivalent to the fact that u = (Aw(t - T))~~ m / 2 e~f solves 
the adjoint heat equation (2.9). 



3 The [RH) a flow as a gradient flow 

In this section, we introduce an energy functional 5" Q for the {RH) a flow, which corresponds 
to Perelman's ^-energy for the Ricci flow introduced in [31, Section 1]. For a detailed study 
of Perelman's functional, we refer to Chow et al. [5, Chapter 5], Miiller [27, Chapter 3], 
or Topping [43, Chapter 6]. In the special case N C R, the corresponding functional was 
introduced by List [22]. We follow his work closely in the first part of this section. 



The energy functional and its first variation 

Let g = gij G T( Sym^_ (T*M)) be a Riemannian metric on a closed manifold M, f : M — >• R 
a smooth function and 4> G C°°(M, N) := {(/) G C°°(M,R d ) | <f>{M) C N}. For a constant 
a(t) = a > 0, we set 

7M,f) ■= J M (R g + \Vf\ 2 g -a\V<f>\iy-fdV g . (3.1) 



Take variations 

gt 3 = 9ij + ehij, hij G l (Sy,n-.:V.\/;). 

f = f + e£, £eC°°(M), 

(p £ = n N (<j) + £1?), ?? G C°°(Af,R d ) with d(x) G T 0(x) iV, 

where 7tat is the smooth nearest-neighbour projection defined on a tubular neighbourhood 
of N C R d . Note that we used the identification T P N C T p R d = R d . We denote by 8 the 
derivative <5 = gj| E _ ) i-e. we have Sg = h, Sf — £ and 5<j> = (dn^ ° (f>)i9 = i9. Our goal is to 
compute 

S3 : a ,g,4,,f(h,0,£) := 4- ? a (g + £h,-K N (<t> + £#),f + ££) 

OLE e=0 

= 8 1 (R+\Vf\ 2 )e-fdV-a-8 f \\7 cf>\ 2 e~ f dV . 

JM JM 

" v ' V v ' 

= : J =: II 

For the variation of the first integral, we know from Ricci flow theory that 



A I 



( /<"(//.., + V<V 3 -/) + (I tr s h-l) (2A/ - |V/| 2 + R))e~ f dV, 
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see Perelman [31, Section 1], or Miillcr [27, Lemma 3.3], for a detailed proof. For the 
variation of the second integral, we compute with a partial integration 

11= [ 2g^\7 l <j) X \7 j i) x e- f dV + [ ( - V^V^ + \V</>\ 2 {\ tr g h - £)) e~ f dV, 
(-2tf A (A g A - (V0 A , V/> fl ) - W ; V# A V,0 A + | (i tr s h-£)) e~UV. 



1 M 

Hence, by putting everything together, we find 



5? a . g ^j(h,$,£) = [ //' ••(/.',, • V,V.,/ - aV^V^e-Uv 

JM 

+ [ (|tr g h - £)(2A f -|V/| 2 + R- a\V(f>\ 2 )e- f dV (3.2) 
Jm 

+ [ 2oe&(T g( j>- {Vcj),Vf})e- f dV, 
Jm 

since ?9A g = $T g (f>, where T g (j> := A g — ^4(0)(V</), V0)m denotes the tension field of </>. 
Gradient flow for fixed background measure 

Now, we fix the measure dfi = dV , i.e. let / = — log (^), where ^ denotes the Radon- 
Nikodym differential of measures. Then, from = Sdfi = (i tr 9 h — t) dfi we get I = | tr g /i. 
Thus, for a fixed measure the functional 5F a and its variation b'S a depend only on g and <f> 
and their variations Sg = h and <5(/> = i9. In the following we write 

S2(ff,0) := J Q ( 5 ,0,-log(^)) (3.3) 

and 

Equation (3.2) reduces to 

S Kg,4 h > = / ( / '"( /,> '.- + V * V j/ - aVi^Vj^) + 2c«?(T g - (V0, V/) (3.4) 

Let {g,4>) € T(Sym 2 + (T*M)) x C°°(M,N) and define on i? := iJ g , = F(Sym 2 (T*A/)) x 
T r f,C 00 (M 1 N) an inner product depending on a and the measure by 

((k ij ,iP),(h ij ,d)) Ha := / (i^%+2a^)d/i. 

From <^ iS> (M) = (gradJ£( 3 ,0), (M)>jr >ai/ , we then deduce 

grades, <£)= (-2(% + ViV,/-aV^V^), r g - (V0, V/) ). (3.5) 

Let 7Ti, 7T2 denote the natural projections of iJ onto its first and second factors, respectively. 
Then, the gradient flow of SFg is 

f =*"i(gradJ£ ( ff , </>)), 
1 ^0 = 7r 2 (grad^( 5 ,0)). 
Thus, recalling the equation J^/ = £ = | tr g ( J^ij), we obtain the gradient flow system 

#.</,.. = ~2(i?y + V,V ,/ - 
&<t> = T g <f>-(V<j>,Vf), (3.6) 
f/ = -i?-A/ + a|V0| 2 . 
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Pulling back with diffeomorphisms 



As one can do for the Ricci flow (see Perelman [31, Section 1] or Miiller [27, page 52]), 
we now pull back a solution (g, 0, /) of (3.6) with a family of diffeomorphisms generated 
by X — V/. Indeed, recalling the formulas for the Lie derivatives (£v/<?)y = 2VjVj/, 
£v/</> = (V</>, V/) and £v// = |V/| 2 , we can rewrite (3.6) in the form 
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2Rc + 2aS7cf> ® V<j> - {L wf g), 

m ( P = T 9^- ( £ v/0), 

I |/ = -A/ + |V/| 2 - R+ a\V</>\ 2 - (£ V //). 

Hence, if ipt is the one-parameter family of diffeomorphisms induced by the vector field 
X(t) = V/(i) e r(TM), i e [0,T), i.e. if ^Vt = o i/> t , V'o = id, then the pulled-back 
quantities g = ip^g, </> = ^<f>, f = ipff satisfy 



JLf 

dtJ 



-2Rc + 2aV0 <g> V0, 



9YJ 

A,/ 



|V/|?-ii + a|V< 



Here, Rc and R denote the Ricci and scalar curvature of g and A, t and the norms are 
also computed with respect to g. In the following, we will usually consider the pulled-back 
gradient flow system and therefore drop the tildes for convenience of notation. 



Note that the formal adjoint of the heat operator □ 



= — t^t — A — T;ti-gh. Indeed, for functions v,w: M x [0, T] 
computation yields 



Jj — A under the flow -j^g = h is 
R, a straightforward 



(Dv)w dVdt 



JM 



dV 



v(D*w) dVdt. 



In our case where h; 



-2R ij + 2aV i (f>Vj<f), this is □* 



M 



- R— a|V '4>\ 2 and thus the 



evolution equation for / is equivalent to e •* solving the adjoint heat equation d*e * = 0. 
The system now reads 

-2Rc + 2aV(f> ® V0, 

9 0, (3.7) 



= □* e - / . 

This means that (RH) a can be interpreted as the gradient flow of 3^ for any fixed back- 
ground measure /i. Moreover, using (3.4), (3.5) and the diffcomorphism invariance of 3q, 
we get the following. 



Proposition 3.1 

Let (g(t), <f>{t))te[0,T) 06 a solution of the {RH) a flow with coupling constant ct(t) = a > 
and let e~^ solve the adjoint heat equation under this flow. Then the energy functional 
SFq, (</, 0, /) defined in (3.1) is non- decreasing with 

^7? a = I (2\Rc-aV<j)®V<j) + licss(f)\ 2 +2a\T g (t)-(\7(t>,X7f)\ 2 )e- f dV>0. (3.8) 
dt J M V / 

Moreover, J a is constant if and only if (g(t), <fi(t)) is a steady soliton. 
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Allowing also time-dependent coupling constants a{t), we obtain the following. 
Corollary 3.2 

Let (g(t), 4>(t))t£[o,T) solve (RH) a for a positive coupling function a(t) and let e~f solve the 
adjoint heat equation under this flow. Then S" a (t) (g(t), (b(t), f(t)) satisfies 



d_ 

dt 



J a = [ (2 |Rc - aV(j) ® + Hess(/)| 2 + 2a \r g <j) - (V<£, V/) | 2 - a\ V<j)\ 2 ) e~ f dV, 

J M ^ ' 

in particular, it is non- decreasing if a(t) is a non-increasing function. 

Minimizing over all probability measures 

Following Perclman [31], we define 

\ a {g^) :=inf {3*(jg,4>) \ ft(M) = 1} = inf {^(s, <t>, /) j e-*W = l^. (3.9) 

The first task is to show that the infimum is always achieved. Indeed, if we set v = e~^/ 2 , 
we can write the energy as 

2 a {jgA,v) = J ( Rv2 + 4|Vtf - a|V0|V)dV = J v(Rv - 4Au - a\\7 <j)\ 2 v^dV. 



Hence 

K(g,<P) = inf 



|y v(Rv-4Av-a\V<j>\ 2 v S )dV J v 2 dV = l 



is the smallest eigenvalue of the operator — 4A + R — a\V<t>\ 2 and v is a corresponding 
normalized eigenvector. Since the operator (for any time t and map <fi(t)) is a Schrodingcr 
operator, there exists a unique positive and normalized eigenvector v m i n (t), see for example 
Reed and Simon [35] and Rothaus [36] . From eigenvalue perturbation theory, we see that if 
g(t) and 4>{t) depend smoothly on t, then so do X a (g(t),(f>(t)) and v m in(t). 

Proposition 3.3 

Let (g(t) , </>(i))te[o,T) be a smooth solution of the (RH) a flow with constant a(t) = a > 0. 
Then \ a (g,<j)) as defined in (3.9) is monotone non- decreasing in time and it is constant if 
and only if 

Jo = Rc-aV0®V0 + Hess(/), 
\o = t 9 0-<V&V/>, 
for the minimizing function f = — 21ogi> m i„. 

Proof. Pick arbitrary times t\, t 2 G [0, T) and let v min (t2) be the unique positive minimizer 
for X a (g(t2),4>(t2))- Put ufa) = v min(^) > and solve the adjoint heat equation D*u = 
backwards on [ti, t^]- Note that u(x' , t') > for all x' G M and t' G [ii, £2] by the maximum 
principle and the constraint J M u dV = J M v 2 dV = 1 is preserved since 

if udV =f (mu)dV+[ u{§- t dV) = -f AudV = 0. 

al JM JM JM JM 
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Here, we used §- t dV = -\tT g {§- t g)dV = {-R + a\V<t>\ 2 )dV and f u = (-A + R - a\V<j)\ 2 )u : 
the latter following from D*u = 0. Thus, with u{i) = e~^*' for all t e [ii , ^2] , we obtain 
with Proposition 3.1 

A„(ff(fi), 4>(h)) < ? a (g(h), <£(ti), /(t x )) < J Q ( 5 (t 2 ), 0(t 2 ), /(t 2 )) = A a (s(t 2 ), 0(t 2 )). 
The condition (3.10) in the equality case follows directly from (3.8). □ 

Again the monotonicity of X a (g, 4 1 ) is preserved if we allow a positive non-increasing coupling 
function a(t) instead of a time-independent positive constant a. 

4 Short-time existence and evolution equations 

Due to diffeomorphism invariance, our flow is only weakly parabolic. In fact, the principal 
symbol for the first equation is the same as for the Ricci flow since the additional term is 
of lower order. Thus, one cannot directly apply the standard parabolic existence theory 
Fortunately, shortly after Hamilton's first proof of short-time existence for the Ricci flow 
in [14] which was based on the Nash-Moser implicit function theorem, DeTurck [10] found 
a substantially simpler proof which can easily be modified to get an existence proof for 
our system (RH) a . Note that we only consider the case where M is closed, but following 
Shi's short-time existence proof in [38] for the Ricci flow on complete noncompact mani- 
folds, one can also prove more general short-time existence results for our flow. We first 
recall some results for the Ricci flow, following the presentation of Hamilton [15] very closely 

Since some results strongly depend on the curvature of (TV, 7) , it is more convenient to work 
with <j> : M — > N itself instead of 4> '■ M — > M. d as in the last section. Therefore, repeated 
Greek indices are summed over from 1 to n = dim N in this section. 

Dual Ricci-Harmonic and Ricci-DeTurck flow 

Let g(t) be a solution of the Ricci flow and : (M,g) — > (M, h) a one parameter family 
of smooth maps satisfying the harmonic map flow JjV = T g x l } with respect to the evolving 
metric g. Note that this is (RH) a =o. If ip(t) is a diffeomorphism at time t = 0, it will stay 
a diffeomorphism for at least a short time. Now, we consider the push-forward g := tp*g of 
the metric g under ip. The evolution equation for g reads 

= -2Ri S + (jc v gh = -2-Rij + ViVj + Vj-y,, (4.1) 

where V denotes the Levi-Civita connection of g and J^ip — — V o if). One calls this the dual 
Ricci-Harmonic flow or also the h-Row. An easy computation (see DeTurck [10] or Chow 
and Knopf [4, Chapter 3]) shows that V is given by 

^ = ff y (*rf 3 .- h r&), (4.2) 

the trace of the tensor which is the difference between the Christoffel symbols of the con- 
nections of g and of h, respectively. Note that the evolution equation of g involves only the 
metrics g and h and not the metric g, and since it involves "r for the fixed background 
metric h it is no longer diffeomorphism invariant. Indeed, one can show (see Hamilton [15, 
Section 6]) that 

&9ij=9 kt VkVm, (4.3) 
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where V denotes the connection of the background metric h. Since V is independent of g 
and V only involves hrst derivatives of g this is a quasilinear equation. Its principal symbol 
is ct(£) = g lj • id, where id is the identity on tensors g. Hence this flow equation is 
strictly parabolic and we get short-time existence from the standard parabolic theory for 
quasilinear equations, see e.g. [24] for a recent and detailed proof. If we additionally assume 
that (M, h) = (M, go) and 0(0) = idM, the flow g which has the same initial data g(0) = go 
as g is called the Ricci-DeTurck flow [10]. 

Now, one can find a solution to the Ricci flow with smooth initial metric <?(0) = go as follows. 
Chose any diffeomorphism 0(0) : M — > M. Since g(Q) is smooth, its push-forward <?(0) is 
also smooth and equation (4.3) has a smooth solution for a short time. Next, one computes 
the vector held V with (4.2) and solves the ODE system 

§- t i> = -Vo^, 

One then recovers g as the pull-back g = ip*9- This method also proves uniqueness of 
the Ricci flow. Indeed, let g l (t) and g 2 (t) be two solutions of the Ricci flow equation 
for t G [0,T) satisfying ^(0) = g 2 (0). Then one can solve the harmonic map heat flows 
= i € {!> 2 i with = V' 2 (0)- Thi s yields two solutions g l = ip\g l of the dual 

Ricci-Harmonic map flow with the same initial values, hence they must agree. Then the 
corresponding vector fields V 1 agree and the two ODE systems -§[tp l = —V l oifj l with the 
same initial data must have the same solutions ip 1 = ip 2 . Hence also the pull-back metrics 
g 1 and g 2 must agree for all t <G [0, T). 

For the dual Ricci-Harmonic flow, the evolution equations in coordinate form, using only 
the fixed connection V of the background metric h, are (see e.g. Simon [39]) 

m9H = g ke V k V e g l3 - ~g M g ip h pq R jkql - g ke g 3P h pq R tkqe ^ 
+ \g M g pq ! (V 'iffpfcV ' 3 g q i + 2Vk9ipV q gjt - 2V k g ip Vtgj q - 4V,g pfe V 'ig jq ), 

where l\, = (Rm(/i))qti denotes the Ricmannian curvature tensor of h. 

Recent work of Isenberg, Guenther and Knopf [17], Schniircr, Schulze and Simon [37] and 
others shows that DcTurck's trick is not only useful to prove short-time existence for the 
Ricci flow, but is also useful for convergence and stability results. Their results show that 
the /(-flow itself is also interesting to study. In this article however, we only use it as a 
technical tool. 



Short-time existence and uniqueness for (RH) a 

Let (g(t), </>(£) )te[o,T) be a solution of the (RH) a flow with initial data (g(0), 0(0)) = (go, 4>o)- 
As for the Ricci-DeTurck flow above, we now let il)(t): (M,g(t)) — > (M,go) be a solution 
of the harmonic map heat flow Jj?/> = r g ^p with ip(0) = id^/ and denote by (g(t),(j>(t)) the 
push-forward of (g(t),cf>(t)) with ip. Analogous to formula (4.1) above, we find 

f = Mm*) + L ^ = T ~90 + <V0, V), 

where V = — 9o r| ■) and V denotes the covariant derivative with respect to g. Note 

that Vi0Vj0 = (d<fi<3d<j))ij as well as (V0, V) = d(f>(V) are independent of the choice of the 
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metric. Using (4.3), we find 

mm = ~9 M ^k^fg lJ + 2a(V4> 8> V0)« , (4.6) 

which is again quasilincar strictly parabolic. The explicit evolution equation involving only 
the fixed Levi-Civita connection V of go can be found from (4.4) by adding 2aVi<^>V J -0 on 
the right and replacing h by go. 

The evolution equation for <j)(t) in terms of V can be computed as follows. Using normal 
coordinates on (N, 7) we have = at the base point and thus T g 4> = A g (j). We find 

§- t 4> x = A^ A + (V0 A , V) = g M (d k d e 4> x - ^V#) + V^V 

= g k \v^^ x + ga n^^ x ~ s n^^ x ) + v# ■ - 9 °n e ) (4.7) 

= g M V k V t 4> X . 

Putting these results together, we have proved the following. 
Proposition 4.1 

Let (g(t), <j)(t)) be a solution of the (RH) a flow with initial data (g(O),0(O)) = (gcb^o)- Let 
ijj{t) : (M,g(t)) — > (M,go) solve the harmonic map heat flow -^ip = T g ijj with tJj(0) = id^/ and 
let (g(t), ^(t)) denote the push-forward of(g(t), 4>(t)) with ?p. Let V be the (fixed) Levi-Civita 
connection with respect to go- Then the dual (RH) a flow (g(t), (f>(t)) satisfies 

§- t 9ij = g k ^k^eg l3 - g u gi P g P o q R]kqi - g M gj P g P o q Rik q i> + 2aV i 4>v j 4> 

+ \g M g pq (yi9pk$ 'jg q i + 2S/ k g w V q gjt - 2V k g ip Veg ]q - ^i7ig pk Vig ]q 



f 4> x = g M V k V^ x + ~g U ( N T x v o 0)V fc ^V/> 



In particular, the principal symbol for both equations is <r(£) = g lj ■ id, i.e. the push- 
forward flow is a solution to a system of strictly parabolic equations. 

Short-time existence and uniqueness for the dual flow (and hence also for the (RH) a flow 
itself) now follow exactly as in the simpler case of the Ricci and the dual Ricci-Harmonic 
flow described above. 

Evolution equations for R, Rc, |V0| 2 and V0 (g> V0 

In the following, we often use commutator identities on bundles like T*M ® (jfTN . The 
necessary formulas arc collected in the appendix. We denote the Riemannian curvature 
tensor on (N, 7) by N Rm and let Rm, Rc and R be the Riemannian, Ricci and scalar 
curvature on (M, g). Moreover, we write 

(Rc, V(f) ® V</>) := RijVi^Vj^, 
(^RmtV^V^V^V^) := X^V^V^V^V^. 

Finally, we use the fact that T g (j> = V p V p </> for the covariant derivative V on T*M (g> <fi*TN, 
cf. Jost [18, Section 8.1]. From the commutator identities in the appendix, we immediately 
obtain for </> € C°°(M,iV) 



+ RipWp^Wjct) + R jp S7 p <f>Vi<f> - 2 ( Ar Rm(V l 0, V p 0)V p 0, V ' 
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Remark. Taking the trace and using T g <j) = 0, we find the well-known Bochncr identity for 
harmonic maps, cf. Jost [18, Section 8.7], 

-A 3 |V<?i| 2 + 2|V 2 0| 2 + 2(Rc,V0® V0) = 2( w Rm(V. t 0, Vj0)V^, V 4 0) . 

Now, we compute the evolution equations for the scalar and Ricci curvature on M. 

Proposition 4.2 

Let (g(t),<p(t)) be a solution to the (RH) a flow equation. Then the scalar curvature evolves 
according to 

§- t R = AR + 2|Rc| 2 - 4a(Rc, V</> <g> V0) + 2a\ Tg (j}\ 2 - 2a|V 2 0| 2 
+ 2a (^RmOv 7 ^, V^)V^, V^) 



and the Ricci curvature evolves by 

>R„w„/?„„ 4- 2rvr„rAV,V,rA - 2rtV„V.-fl l )V„ V a6 

(4.10) 



§- t Rij = A L R l3 - 2R iq R jq + 2R ipjq R pq + 2ar g (j)X7 i X7 j (j) - 2aV p V^VpV j( 
+ 2aRpi jq V p <l>V q <f> + 2a ( N Rm{\/ l( j), V p 0)V p 0, Vrf) . 



Here, Al denotes the Lichnerowicz Laplacian, introduced in [21], which is defined on sym- 
metric two-tensors Uj by 

Aljtij . — Atij -(- 2Ripj q tp q Riptpj Rjptpi- 

Proof. We know that for Jj<7ij — hij the evolution equation for the Ricci tensor is given by 

f t Rij = -\A L hij + \V % Vph PJ + \VjVphpi - |ViV,-(tr fl h), (4.11) 

see for example [27, Proposition 1.4] for a proof of this general variation formula. For 
hij = —2Rij, we obtain (with the twice contracted second Bianchi identity) j^Rij = A^i?^-. 
For h^ = 2V,0Vj0, we compute, using (4.8) 

§- t Rij = -A L {Vi4>Vj4>) + ViV p (Vp^Vj0) + V J V p (Vp0V l 0) - ViV^V^Vj,^) 

= 2T g 4>V i V j <t> - 2V i Vp4>V j Vp<j) - 2R vpjq V P 4>V q (j> + 2 ( ff Rm(V^, Vp0)Vp</>, Vrf) ■ 

Linearity then yields (4.10). For the evolution equation for R, use 

§- t R = ^Rn) = + (2% - 2aV l 0V J 0)%. 

The desired evolution equation (4.9) follows. An alternative and more detailed proof can be 
found in the author's thesis [28, Proposition 2.3]. □ 

Next, we compute the evolution equations for |V^| 2 and V</> (g> V</>. 
Proposition 4.3 

Let (g(t), (j)(t)) be a solution of (RH) a . Then the energy density of (j) satisfies the evolution 
equation 

£ t \V<p\ 2 = A|V(/>| 2 - 2a|V0 g> V0| 2 - 2|V 2 ^| 2 + 2 (^RmfV^, VjflVrf, Vi<f>) . (4.12) 
Furthermore, we have 



f t (V^V^) = A(V^V^) - 2VpVi^VpV^ - RipVp4>S7 j( j) - R JP Vp4>Vi 
+ 2 ( N Rm(V^, V p 0)Vp0, V j( />) . 



(4.13) 
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Proof. We start with the second statement. We have 

i(V^V^) = (VtV^)V^ + (VtV^)V^ 
= V < (^)V J ^ + V J -(^)V i 

= ViTg<j)Vj<f> + V jTgffi i<j>, 

where the meaning of the covariant time derivative V t is explained in the appendix. The 
desired evolution equation (4.13) now follows directly from (4.8). We obtain (4.12) from 
(4.13) by taking the trace, 

f |V0| 2 = (2R ij - 2aV^V^)V^V^ + 3 ij |(V^V^) 

= -2a|V0® V0| 2 + A|V</f - 2|V 2 0| 2 +2( Ar Rm(V^,V^)V 3 </>,V l 0) . □ 



Evolution of S = Rc — aV</> £g> V</> and its trace 

We write again § := Rc — aS/cj) ® V(f> with components Sij = Rij — aVj^V j(f> and let 
S — R — a\V<t>\ 2 be its trace. Then, we can write the (RH) a flow as 

{~§i9ij = 
m<i> = T <^ 

and the energy from Section 3 as 7 a {g, 0, /) := J M (5+ \\7 f\ 2 )e~? dV g . It is thus convenient 
to study the evolution equations for S and S. Indeed, many terms cancel and we get much 
nicer equations than in the previous subsection. 

Theorem 4.4 

Let (g(t),4>(t)) solve (RH) a with a(t) = a > 0. Then S and S defined as above satisfy the 
following evolution equations 

£S = AS + 2|^| 2 + 2a|r^| 2 , 
feSy =A L S ij + 2ar g <jN i W j <t>. 

Proof. This follows directly by combining the evolution equations from Proposition 4.2 with 
those from Proposition 4.3. □ 

Remark. Note that in contrast to the evolution of Rc, R, £g> V</> and |V</>| 2 the evolution 
equations in Theorem 4.4 for the combinations Rc — aV</> <8> Vcj> and R — a|V</>| 2 do not 
depend on the intrinsic curvature of N. 

Corollary 4.5 

For a solution (g(t), 4>(t)) of (RH) a with a time- dependent coupling function a(t), we get 

§- t S = AS + 2\S. tJ \ 2 + 2«M| 2 - d|V</>| 2 , 
■jftSij = A L Sij + 2a.T g 4>ViV j<p - 6tWi<t>Wj4>. 

5 First results about singularities 

In this section, we often use the weak maximum principle which states that for parabolic 
partial differential equations with a reaction term a solution of the corresponding ODE yields 
pointwisc bounds for the solutions of the PDE. Since we work on an evolving manifold, we 
need a slightly generalized version. The following result is proved in [4, Theorem 4.4]. 
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Proposition 5.1 

Let ii:Mx[0,T)-)Ife(i smooth function satisfying 

§- t u > A g(t) u + (X(t),Vu) g(t) + F(u), (5.1) 

where g(t) is a smooth 1-parameter family of metrics on M , X(t) a smooth 1-parameter 
family of vector fields on M , and F : R — > R is a locally Lipschitz function. Suppose that 
m(-,0) is bounded below by a constant Cq € R and let <f>(t) be a solution to 

&4> = F(<P), 0(0) = C . 

Then u(x,t) > <j)(t) for all x £ M and all t G [0,T) for which (j)(t) exists. 



Similarly, if (5.1) is replaced by j^u < A g r t \u + (X(t), Vu) g u\ + F(u) and u(-, 0) is bounded 
from above by Co, then u(x,t) < <fi(t) for all x G M and t € [0,T) for which the solution 
4>(t) of the corresponding ODE exists. 



Using this, an immediate consequence of Corollary 4.5 is the following. 



Corollary 5.2 

Let (g(t), 4>(t)) be a solution to the (RH) a flow with a nonnegative, non-increasing coupling 
function a(t). Let S(t) = R(g(t)) — a(t)\V 4>{t)\ 2 g ^ as above, with initial data 5(0) > on 
A/. Then R m in(t) ■= min xe M R(x, t) — > oo in finite time and thus g(t) must become singular 
in finite time T S i ng < „„ m < oo. 



Proof. Since a(t) > and a(t) < for all t > 0, Corollary 4.5 yields 

f S> AS + 2\Sij\ 2 > AS+^S 2 (5.2) 
and thus by comparing with solutions of the ODE -^a(t) = ^a{t) 2 which are 

q(0) 

°W = i_as o (0)' 

the maximum principle above, yields 

5mm(t) - 1 -if2(o) (5 - 3) 

for all £ > as long as the flow exists. In particular, if S m i n (0) > this implies that 
S m in(t) — ^ oo in finite time To < 2S m < oo. Since R — S + a|V0| 2 > S, we find 
that also R m i n (t) — > oo before To and thus g(t) has to become singular in finite time 

T smg <T < 25 m 7„(0) < 00 ■ ^ 

As a second consequence, we see that if the energy density e(<p) = ^|V(/)| 2 blows up at 
some point in space-time while a(t) is bounded away from zero, then also g(t) must become 
singular at this point. 



Corollary 5.3 

Let (g(t), <fi(t))te[o,T) be a smooth solution of (RH) a with a non-increasing coupling function 
a(t) satisfying a(t) > a > for all t G [0, T). Suppose that \\74>\ 2 (xk, tk) — >• oo for a sequence 
(xk,tk)keN with tk T. Then also R(xk,tk) — > oo for this sequence and thus g(tk) must 
become singular as tk approaches T. 



17 



Proof. From (5.2) we obtain 5 > S mm (0) and thus 

# = a|V0| 2 + 5>g|V0| 2 + 5 mm (O), V(x, t) G M x [0, T). (5.4) 

Hence, if \\7(f>\ 2 (xk, tk) — > oo for a sequence (xk,tk)keN C M x [0, T) with ^ /" T then also 
R(xk,tk) — > oo for this sequence and g(tfe) must become singular as t;t /T. □ 

Remark. The proof shows that Corollary 5.3 stays true if a(t) \ as i / T as long as 
\V4>\ 2 (xk,tk) — > oo fast enough such that a(tk)\ V(f>\ 2 (xk, tk) — > oo still holds true. 

Now, we derive for t > an improved version of (5.4) which does not depend on the initial 
data 5(0). Using (5.2) and the maximum principle, we see that if S m i n (0) > C £ I we 
obtain 

m 

and thus 5(i) > — S for all t > as long as the flow exists, independent of 5(0). More 
rigorously, this is obtained as follows. The inequality (5.2) implies 

£ t (tS) = S + t(£ t S)>A(tS) + S(l + %S). 

If (xq, to) is a point where tS first reaches its minimum over M x [0, T— S], 5 > arbitrarily 
small, we get 5(a;o,*o)(l + ^fS(xo,to)) < 0, which is only possible for t S(xo,t ) > 
Hence > -y on all of M x [0, T — 5]. Since S was arbitrary, we obtain the desired 
inequality 5(t) > — § everywhere on M x (0,T). This yields 

R > a|V0| 2 - g > «|V0| 2 - § , V(x,t) G M x (0,T), 

which immediately implies the following converse of Corollary 5.3. 

Corollary 5.4 

Let (g{t),<p(t))t^[o,T) be a smooth solution of (RH) a with a non-increasing coupling function 
a{t) > a > for all te [0, T). Assume that R < R on M x [0, T). Then 

|V0| 2 <^ + ^-, V(x,t) € M X (0,T). (5.5) 



Singularities of the type as in Corollary 5.3, where the energy density of <f> blows up, can 
not only be ruled out if the curvature of M stays bounded. There is also a way to rule them 
out a-priori. Namely, such singularities cannot form if either N has non-positive sectional 
curvatures or if we choose the coupling constants a(t) large enough such that 

Here K denotes the sectional curvature of N. More precisely, we have the following 
estimates for the energy density |V0| 2 . 

Proposition 5.5 

Let (g{t), </>(t))ie[o,T) be a solution of (RH) a with a non-increasing a(t) > and let the 
sectional curvature of N be bounded above by N K < cq. Then 
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i) if N has non-positive sectional curvatures or more generally if c — < 0, the energy 
density of <p is bounded by its initial data, 

|V#M)| 2 < max|V^(y,0)| 2 , V(x,t) G M x [0,T). (5.6) 

3/SM 

ii) if N has non-positive sectional curvatures and a(t) > a > 0, we have in addition to 
(5.6) the estimate 

|V#M)| 2 <£^, V(z,t)eMx(0,T). (5.7) 

Hi) in general, the energy density satisfies 

\V<j)(x,t)\ 2 < 2max|V</>(?/,0)| 2 , V(x, t) G M x [0, T*), (5.8) 

where T* := minjT, (4c max^^l V^(y, 0)! 2 )" 1 } . 

Proof. This is a consequence of the evolution equation (4.12) and the Cauchy-Schwarz in- 
equality 

i|V^| 4 = £| 5 y Vi0V^| a < |V,0V^| 2 < |V0| 4 . (5.9) 

i) If N has non-positive sectional curvatures, ( Ar Rm(Vi0, VjfyVjcj), Vitp) < 0, the evo- 
lution equation (4.12) implies 

^|V0| 2 < A|V0| 2 . (5.10) 

If co - 2® < 0, we have 2( A fem(V i ^V J 0)V j 0,V i ^) < 2c o |V0| 4 < 2^|V0| 4 < 
2a|Vi</>Vj</>| 2 , and we get again (5.10) from (4.12). The claim now follows from the 
maximum principle applied to (5.10). 

ii) If (^Rm^i^Vj^Vj^V^) < 0, (4.12) and (5.9) imply 

f |V0| 2 < A|V0| 2 - 2a|V 4 0V,0| 2 < A|V0| 2 - 2^|V0| 4 . 

We obtain 

|(t|W>| 2 ) = |V0| 2 +t(f t |V0| 2 ) < A(t|V0| 2 ) + |V0| 2 (1 - 2t^|V0| 2 ). 

At the first point (xo , io ) where i 1 V ^ | 2 reaches its maximum over M X [0 , T — 5] , <5 > 
arbitrary, we find 1 — 2< ^ |V0| 2 (xo, £o) > 0, i.e. 

io|V0| 2 (x o ,<o) < ^ < ^, 

which implies that t|V</>| 2 < g for every (x, t) G M X [0, T - 6]. The claim follows. 

iii) From (4.12), we get 

f |V^| 2 < A|V0| 2 + 2c o |V0| 4 . 
By comparing with solutions of the ODE = 2coa(t) 2 , which are 

a(0) 1 



1 - 2c o a(0)i' ~ 2c o a(0)' 
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the maximum principle from Proposition 5.1 implies 

+W s max y g M |V(/>(y,0)| 2 
|V ^ MI ^l-2c o max yeM |V^0)Pt' ^ 

for all x € M and t < minjT, (2co max yS M| V</>(j/, 0)| 2 ) -1 }. In particular, this proves 
the doubling-time estimate that we claimed. □ 

6 Gradient estimates and long-time existence 

For solutions (g(t), 4>(t)) of the (RH) a flow with non-increasing a(t) > a > 0, we have seen 
in Corollary 5.4 that a uniform bound on the curvature of (M, g(t)) implies a uniform bound 
on |V</>| 2 . Therefore, one expects that a uniform curvature bound suffices to show long-time 
existence for our flow. The proof of this result is the main goal of this section. 

Evolution equations for Rm and V 2 (j) 

With J^ffy = hij := — 2i?y +2aVi0 M Vj0 M , we find the evolution equation for the Christoffcl 
symbols 

= g pq {-V l R 3q - S7 3 R iq + VqRij) + 2aV l V J </>V J > 

With this, an elementary computation yields the following evolution equation for the Rie- 
mannian curvature tensor (see [28, Proposition 3.2] for a detailed proof). 

Proposition 6.1 

Let (g(t),<j)(t))t£[o,T) be a solution of (RH) a . Then the Riemann tensor satisfies 



StRijki = ViVfci?^ — ViVfRjk — VjVkRu + VjVeRik ~ RijqiRkq ~ RijkqReq 
+ 2a(V l V fe 0V j V e (j> - ViV^VjV fe <£ - ( JV Rm(V 4 </>, V^)V fc & V^> ) . 



(6.2) 



Remark. Taking the trace of (6.2), we obtain (4.10), using the twice traced second Bianchi 
identity. This gives an alternative proof of Proposition 4.2. 

If we set a = in (6.2), we obtain the evolution equation for the curvature tensor under the 
Ricci flow. It is well-known that this evolution equation can be written in a nicer form, in 
which its parabolic nature is more apparent. In [14, Lemma 7.2], Hamilton proved 



(6.3) 



ViV^Rjt — ViV (Rjk — VjVkRu + VjVfi?^ 

= ARijke + 2(Bijke — Bijik — B^jk + Bikje) — RpjkeRpi — RipkeRpj, 

where Bijki '■= RipjqRkpeq- Plugging this into (6.2) yields the following corollary. 
Corollary 6.2 

Along the (RH) a flow, the Riemannian curvature tensor evolves by 

-§tRijke = AR^ke + KBijke — B^ik — Bujk + Bikje) 

— (RpjkeRpi + RipkeRpj + Rij P eR P k + RijkpRpe) (6.4) 
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There is a useful convention for writing such equations in a short form. 



Definition 6.3 

For two quantities A and B, we denote by A * B any quantity obtained from A <g> B by 
summation over pairs of matching (Latin and Greek) indices, contractions with the metrics 
g and 7 and their inverses, and multiplication with constants depending only on m = dim M , 
n = dim A and the ranks of A and B. We also write (A)* 1 := 1 * A, (A)* 2 = A* A, etc. 



This notation allows us to write (6.1) and (6.4) in the short forms 

§- t T= (VRm)* 1 +aV 2 0* V0. (6.5) 

and 

|Rm = ARm + (Rm)* 2 + a(V 2 0)* 2 +a N Rm * (V</>)* 4 . (6.6) 

It is now easy to compute the evolution of the length of the Ricmann tensor. Together with 
Jj<? _1 = (Rm)* 1 + a(V0)* 2 , the above formula yields 

§- t |Rm| 2 = {§- t g- 1 ) * Rm * Rm + 2R iM {§ i R ijH ) 

= A|Rm| 2 - 2|VRm| 2 + (Rm)* 3 + «(Rm)* 2 * (V0)* 2 (6.7) 
+ qRie * (V 2 </))* 2 + aRm * N Rm * (V</))* 4 . 

Corollary 6.4 

Along the (RH) a flow, the Riemannian curvature tensor satisfies 



^|Rm| 2 < A|Rm| 2 -2|VRm| 2 + C|Rm| J + aC|V</>| 2 |Rm|' 
+ aC|V 2 </>| 2 |Rm| + aCc o |V0| 4 |Rm|, 



(6.8) 



for constants C > depending only on the dimension of M and Cq = cq(N) > depending 
only on the curvature of N . If N is flat, we can choose Cq = 0. 



Proof. Follows directly from (6.7) and the fact that | w Rm| is bounded on compact A. □ 

For the evolution equation for the Hessian of <j), it is important that we do not use the 
*-notation directly. Indeed, we will see that all the terms containing derivatives of the cur- 
vature of M cancel each other (using the second Bianchi identity), a phenomenon which 
cannot be seen when working with the *-notation. 

A short computation using (A.l) and (A. 2) shows that the commutator [ViVj,A]0 A = 
Vi\7jT g (f> x - AViVj(j> x is given by 

[ViVj, A]0 A = V k R mik V p (j) X + 2R lk]p V k V p <j) X 

- RipVj ;V P 4> X - V ? i? JP V p A - RjpVtVp^ (6.9) 
+ ( w Rm * V 2 ^ * (V0)* 2 + (a w Rm) * (V0)* 4 ) .., 

see [28, equation (3.10)] for details. With (6.1) we continue 

[ViVj, A]</> A - (^r^.)Vfc^ = (Rm * V - 2aV i V,0V^V fc </> A 

+ ( w Rm * V 2 (f> * (V</>)* 2 + (d w Rm) * (V(/))* 4 ) i . , 
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where we used the second Bianchi identity (VkRjpik + VjRip — V p i?y)V p A — to cancel 
all terms containing derivatives of the curvature of (M,g). Since the V 2 </> live in different 
bundles for different times, we work again with the covariant time derivative Vt (and with 
the interpretation of V 2 </> as a 2- linear TN- valued map along <fi), as we already did in Section 
4, see appendix for details. At the base point of coordinates satisfying (A. 5), we find with 
(A. 6) and the remark following it 



V t (V 4 V^ A ) = ViVj^ - (<|lt,)V fc </> A + w Rm(f <f>, V^)Vj 



A 

= AViVj0 A + (Rm * V 2 A )*j + aViVj^ * V0 * V0 A (6.10) 
+ (%! * V 2 * (V0)* 2 + (d^Rm) * (V0)* 4 ) . 

With A|V 2 0| 2 = 2A(V»V i (/) A )V i V i (/) A + 2|V 3 c/>| 2 , we finally obtain 

f |V 2 0| 2 = (f t g- 1 ) * (V 2 <£)* 2 +2V t (V l V^ A )V. t V,0 A 

= Rm * (V 2 0)* 2 + a(V0)* 2 * (V 2 0)* 2 + A|V 2 0| 2 - 2|V 3 0| 2 (6.11) 
+ N Rm * (V 2 0)* 2 * (V0)* 2 + (d N Rm) * (V 2 0) * (V0)* 4 

Since | ^Rm) and Ic^Rml are bounded on compact manifolds N, say by a constant c\, this 
proves the following proposition. 

Proposition 6.5 

Let (g{t), 0(i))te[o,T) 6e a solution of {RH) a . Then the norm of the Hessian of <f> satisfies 
the estimate 

| |V 2 0| 2 < A|V 2 ^| 2 - 2|V 3 0| 2 + C|Rm||V 2 0| 2 

+ aC|V</>| 2 |V 2 </>| 2 + Cci|V0| 4 |V 2 0| + Cci|V0| 2 |V 2 0| 2 

along the flow for some constants C = C(m) > and ci = Ci(N) > depending on the 
dimension m of M and the curvature of N , respectively. If N is flat, we may choose C\ = 0. 

Remark. If we set a = 0, Corollary 6.4 and Proposition 6.5 yield the formulas for the Ricci- 
DcTurck flow (RH)q, in particular (6.8) reduces to the well-known evolution inequality 

^|Rm| 2 < A|Rm| 2 - 2|VRm| 2 + C|Rm| 3 

for the Ricci flow, first derived by Hamilton [14, Corollary 13.3]. Moreover, if a = 2 and 
TV C R (and thus cq = c\ = 0), the estimates (6.8) and (6.12) reduce to the estimates found 
by List (cf. [22, Lemma 2.15 and 2.16]). 

Interior-in-time higher order gradient estimates 

Using the evolution equations for the curvature tensor and the Hessian of </>, we get evolution 
equations for higher order derivatives by induction. 

Definition 6.6 

To keep the notation short, we define for k > 



I k := V i Rm * V J Rm + a ^(<9 l N Rm + 1) * V n <f> * . . . * V je 

i+j=k A k 

+ a^2\/ n cj) * . . . * V^- 1 cj) * V jl Rm, 



B, 



22 



where the last two sums are taken over all elements of the index sets defined by 

A k := \ 0<i<k + l,l<j s <k + 2Vsandj 1 + ... + j e = k + 4}, 

B k := {(ji, ...,j e ) | 1 < j s < k + 2 Vs < 4 < j t < k and jx + . . . + u = k + 2}. 

Lemma 6.7 

Let (g(t), 0(i))ie[o,T) ^ e a so ^ u ti° n to the (RH) a flow. Then for k > ; m'tt defined as in 
(6.13), we obtain 

fV^Rm = AV fc Rm + I k . (6.14) 



Remark. If (iV, 7) = (R, 5) , all terms in I k containing d l %m vanish, and the result reduces 
to (a slightly weaker version of) List's result [22, Lemma 2.19]. Note that we do not need 
all the elements of Ak here, but defining Ak this way allows us to use the same index set 
again in Definition 6.8. 

Proof. From (6.6), we see that (6.14) holds for k = 0. For the induction step, assume that 
(6.14) holds for some k > and compute 

f t V fc+1 Rm = §- t (<9V fc Rm + V * V fc Rm) = V(AV fc Rm) + VI k + §- t T * V fc Rm. 

Since Vifc is of the form Ik+i and also JjF * V fc Rm = (VRm + aV 2 * V(/>) * V fc Rm appears 
in Ik+i, it remains to compute the very first term. With the commutator rule (A.l), we get 

V(AV fc Rm) = AV fc+1 Rm + VRm * V fe Rm + Rm * V fc+1 Rm 

= AV fc+1 Rm + 4 +1 . □ 

Similar to (6.7), we obtain 

j- t |V fc Rm| 2 = (^T 1 ) * V fc Rm* V fc Rm + 2V fe Rm(f V fc Rm) 

= Rm * (V fc Rm)* 2 + a (V</>)* 2 * (V fc Rm)* 2 + 2V fc Rm(AV fe Rm) + V fc Rm * I k . 

Hence, using the fact that Rm * V fc Rm as well as a (V(/>)* 2 * V fe Rm are already contained 
in Ik, we find 

^|V fc Rm| 2 = A|V fe Rm| 2 - 2|V fc+1 Rm| 2 + V fe Rm * I k . (6.15) 

Definition 6.8 

To compute the higher order derivatives of <f>, we define 

J k := V<Rm * ^ j+2( t ) + ^2( di NRm + 1) * V jl ^ * . . . * V Jl 4) 



i+j=k 



(6.16) 



with A k and B k defined as in Definition 6. 6. 



Lemma 6.9 

Let (g{t), 4>(t)) t £[o t T) be a solution to the (RH) a flow. Then for k > 0, with J k defined as 
in (6.16), we have 

V t (V fc+2 0) = AV fe+2 + J k . (6.17) 
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Proof. For k = 0, the statement holds by (6.10). For the induction step, we use again the 
interpretation of X7 k 4> as a fc-linear TAT- valued map along d> and compute analogously to 
(A. 6) and the remark following it 

v t (v fc+3 0) = vv t (v fc+2 0) + f r * v k+2 4> + N RM-§- t & v0)v fe+2 (/) 

= V(AV fc+2 0) + V J k + §- t T * W k+2 (f> + N Rm * V 2 (f> * V</> * V fe+2 0. 

Again, we only have to look at the first term, since VJ&, "Rm * V 2 * V0 * V fe+2 c/> and 
* V fc+2 </> = (VRm + aS7 2 4>* V4>) * V k+2 4> are obviously of the form J fc+1 . With a higher 
order analog to (A. 3), we obtain 

V(AV fc+2 0) = Vj,VVj,V fe+2 <£ + Rm * V fe+3 9i + w Rm * W k+3 <f> * V</» * V0 
= VpV p V fe+3 + VRm * V fe+2 </> + Rm * V fc+3 

+ (<9 w Rm) * V fe+2 * (V(/.)* 3 + "Rm * V fc+3 * (V?!.)* 2 

+ w Rm * V fc+3 c/> * V 2 </> * V</> 
= AV H3 ^+J H1 , 

and the claim follows. □ 
As in (6.11), we compute 

JL |v fc+2 0| 2 = (f y- 1 ) * V fe+2 * V fe+2 + 2V fe+ Vv t (V fc+ V) 
= Rm * (V fe+2 ^)* 2 + a (V0)* 2 * (V fc+2 0)* 2 

+ 2V fc+2 A (AV fc+2 A ) + V fe+2 </> * J fc 
= 2V fe +V(AV fc+ V) + V fc+2 * Jfe. 

With A\V k+2 (/)\ 2 = 2V fc+2 A (AV fc+2 A ) + 2|V fe+3 0| 2 , we finally find 

j>- t \V k+2 cp\ 2 = A|V fe+2 <£| 2 - 2|V fe+3 </>| 2 + V fe+2 </> * J fc . (6.18) 

The next trick is to combine the two equations (6.15) and (6.18) to a single equation. 
Remember that we already used a similar idea in Section 4, where we combined the evolution 
equations of Re and V (/><£> V</> (respectively i? and |V</>| 2 ) to a single equation for a combined 
quantity 5ij (respectively S), which was much more convenient to deal with. Here, we define 
the "vector" 

7 = (Rm, V 2 0) G T((T*M) m ) x r((T*M)® 2 ® 0*TiV) (6.19) 

with norm |T| 2 = |Rm| 2 + |V 2 ^| 2 and derivatives V fe T = (V fe Rm, V fe+2 0). Combining the 
evolution equations (6.15) and (6.18), we get 

f |V fe T| 2 = A|V fc T| 2 - 2|V fe+1 T| 2 + L k , (6.20) 

where Lk '■= V fc Rm * Ik + V k+2 4> * Jk- We can now apply Bernstein's ideas [2] to obtain 
interior-in-time estimates for all derivatives |V fe T| 2 via an induction argument. For the Ricci 
flow, this was independently done by Bando [1] and Shi [38]. 

Theorem 6.10 

Let (g(t), </>(i))te[o,T) solve (RH) a with non-increasing a(t) G [a, a], < a < a < oc and 
T < oo. Let i/ie Riemannian curvature tensor of M be uniformly bounded along the flow, 
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|Rm| < i?o- Then there exists a constant K = K(a, a, Ro,T,m, N) < oo such that the 
following two estimates hold 

|V</)| 2 <f, V(x,t)eMx(0,T), (6.21) 

|T| 2 = |Rm| 2 + |V 2 (/>| 2 < W(x,t) e MX (0,T). (6.22) 
Moreover, there exist constants Ck depending on k, a, m and N , such that 

|V fc T| 2 = |V fc Rm| 2 + |V fc+2 (/)| 2 <C fc (f) fc+2 , V(x,t) eMx (0,T). (6.23) 

Proof. Since the method of proof is quite standard, we only give a brief sketch of the 
argument and refer to the authors thesis [28, Theorem 3.10] for more details. Setting 
1 < A'x := max { 2 "' 2j ^ T+m | R a T, 1} < oo, we obtain 

|V0| 2 < + £L < and |Rm|<^i, V(x, i) G M x (0, T) 

a 2at t t 

from Corollary 5.4. In the following, C denotes a constant depending on Ki, a, m and the 
geometry of N, possibly changing from line to line. With the estimates for |Rm| and IV^' 2 
and using \V 2 (f>\ < ± + i|V 2 0| 2 , we obtain for f(x,t) := t 2 \X7 2 (j)\ 2 (8K 1 + t\V(, 



\2\ 



(f - A)/ < -2t 2 |V 3 </)| 2 (8A' 1 + t\Vcj>\ 2 ) + f / + f • 9#i - 2t 3 |V 2 </>| 4 + S / 
+ St 3 |V 3 ^»||V 2 0| • |V 2 0||V0| 

on A/ x (0,T). The last term can be absorbed by the two negative terms, 

8i 3 |v 3 </>||v 2 </>| 2 |v</>| < i(8^i)(4t 2 |v 3 0| 2 ) + Kstfir^iet^vVHWl 2 ) 

= 2t 2 |V 3 0| 2 -8K 1 + M5C-t 3 |V 2 0| 4 
< 2t 2 |V 3 </)| 2 (8A'i+t|V0| 2 ) +i 3 |V 2 0| 4 . 

Here, we used 8t g^ < 1 which motivates our choice of the constant 8Ki in the definition 
of /. From (f - A)/ < f / + f - i 3 |V 2 </)| 4 < T §^ ?i (C/ + C - / 2 ), we conclude, 
using /(-,0) = and the maximum principle, that — f 2 + Cf + C > 0. Equivalently, 
/ < D := \ (C + VC 2 +4C) on M x [0,T). For positive t, this implies 



^- ^iwi ^nf' ' (6 - 24) 



where A"2 := \J D/8K\ < oo. Setting if := K1+K2, we get (6.21) and (6.22). Using a similar 

argument, one can then prove (6.23) inductively. The crucial estimates are L k < C(-j-) fe+3 
and 

L k+1 <C(f) k+4 + Cf |V fe+1 T| 2 , 

where C now denotes a constant depending only on a, m, iV and fc (but not on K or T). 
Defining h(a?,t) := i fe+3 |V' £+1 T| 2 (A + t fc+2 |V fc T| 2 ) with A = 8C k K k+2 , these estimates give 



(f - A)/i < -2i fc+3 |V fe+2 T| 2 (A + t fc+2 |V fc T| 2 ) + + f A fe+4 (A + i fc+2 |V fc T| 2 ) 
- 2t 2fc+5 |V fc+1 T| 4 + Sa^+V^IV^TI 2 

+ 8 ^+5| V fc+2 T || V fc+l T | . | V fc + i T || V fc T ^ 
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Using K > 1, the inductive assumption and Cauchy-Schwarz, we rewrite this as 
[O. - A)h < -2i fc+3 |V fc+2 T| 2 (A + i fc+2 |V fc T| 2 ) - §i 2fc+5 |V fc+1 T| 4 

Again, the last term can be absorbed by the negative terms 

8< 2fc+5| V fc+2 T || V fc+i T | 2 | V fc T | < i A (4t fc+3 |V fe+2 T| 2 ) + |A- 1 (l6t 3fe+7 |V fe+1 T| 4 |V fe T| 2 ) 

= 2< fc+3 |V fc+2 T| 2 ■ A + Stk+2 ^\ 2 . t 2fc+5| V fc+L T |4 

< 2t fc+3 |V fc+2 T| 2 (A + t fc+2 |V fe T| 2 ) + i 2fc+5 |V fc+1 T| 4 , 

which explains our choice of A. A maximum principle argument like the one for / above 
then yields -h 2 + CK 2k+5 h + CK ik+w > 0, i.e. h < \{C + VC 2 + AC)K 2k+b =: DK 2k+5 
on M x [0,T). For t > 0, 

fc+i ,2 _ h DK 2k+5 (K\k+3 

|V T| ~ t k + 3 {\ + t k + 2 \V k 7\ 2 ) ~ t k + 3 8C k K k + 2 ~ k+1 \~) ' 

where Ck+i ■= D/(8Ck)- This proves the induction step and hence the theorem. □ 

In the following corollary, we state a local version of the gradient estimates. The setting is 
made in such a way to perfectly fit the proof of the non-collapsing result in Section 8. 

Corollary 6.11 

Let {g(t), </>(i))te[o,r) solve (RH) a with non-increasing a(t) G [a, a], < a < a < oo and 
T' < T < oo. Let B := B g (<j<<){x,r) be a ball around x with radius r, measured with respect 
to the metric at time T' . Assume that |Rm| < Rq on the set B x [0,T"). Then there exist 
constants K = AT(g, a, i?o, T, m, N) < oo and Ck = Cf~(k, a, m, N) for k G N, Co = 1, such 
that the following estimates hold for k > 

|V0| 2 < f and |Rm| < f , V(x, t) £ B 1 ^ 2 x (0, T'), (6.25) 
|V fc T| 2 = |V fe Rm| 2 + |V fe+2 0| 2 <C fe (f) fe+2 , V(x,t) G B 1 / 2 x (0,T'), (6.26) 

where B 1 ! 2 := B g ^ T ij(x,r /2) is £/ie &aH o/ /lay i/ie radius and the same center as B. 

Proof. The statement (6.25) follows exactly as in Theorem 6.10. The induction step is 
carried out using a cut-off function to ensure that the maxima are attained in the interior 
of the set B. More details can be found in the authors thesis [28, Proposition A. 5]. □ 

Long-time existence 

This subsection follows Section 6.7 about long-time existence for the Ricci flow from Chow 
and Knopf's book [4] . We first need a technical lemma. 

Lemma 6.12 

Let {g(t),<f>(t)) t £\o t T) solve (RH) a with a non-increasing a(t) G [a, a], < a < a < oo and 
T < oo. Let the Riemannian curvature tensor of M be uniformly bounded along the flow, 
|Rm| < R , and fix a background metric g. Then for each k > there exists a constant Ck 
depending on k, m, N, T, a, a, Rq and the initial data (<?(0), (f>(0)) such that 

\V k g(xM + \f k Rm(xM + \V k ^(xM<C k (6.27) 
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for all (x,t) £ M x [0,T). Here, V = ff V denotes the Levi-Civita connection with respect to 
the background metric g. 

Proof. With Theorem 6.9, the proof becomes a straight forward computation, and we there- 
fore only give a sketch. In [4, Section 6.7], all the details are carried out in the case of the 
Ricci flow and they can easily be adopted to our flow. Since M is closed, there exists a finite 
atlas for which we have uniform bounds on the derivatives of the local charts. Working in 
such a chart ip : U — > M. m , it suffices to derive the desired estimates for the Euclidean metric 
(5 in U and the ordinary derivatives, since g and V are fixed. In particular, we can interpret 
r as a tensor, namely r = T — S T. On the compact interval [0,T/2], all the derivatives 
|V fc </>|g and |V fe Rm|g are uniformly bounded. On the interval [T/2,T), Theorem 6.10 above 
yields uniform bounds for these derivatives. Hence 

IV^Ig + |V fe Rm|2 < C k (6.28) 

for some C k < oo. In particular, S = Rc — aVcj) <£) V0 is uniformly bounded on [0, T). From 
[4, Lemma 6.49], we infer that all g(t) are uniformly equivalent on [0,T), and thus for some 
constant C 

C^S < g(x, t) < C5, V(.t, t) eU x [0, T). (6.29) 
With §- t {dg) = d{§- t g) = -2d§ = -2(V§ + T * S), we compute 

\gdg\s < C\§- t dg\ < C|VS| + C|r||S|. (6.30) 

Thus, we have 

|^r| < C|VRc|+2a|V0||V 2 0| 

which yields a bound for |T| by integration. Together with the bounds for |S| and |VS| that 
we obtain from (6.28), we conclude that \-^dg\s is uniformly bounded, and hence - again by 
integration - \dg\$ is uniformly bounded on U x [0, T). Finally, using a partition of unity for 
the chosen atlas, we obtain a uniform bound for \Vg\g on M x [0,T). A short computation 
as for (6.30) yields 

fc fc-i 
liV% < C\§- t V k g\ < £ Ci |lf IV^'SI + J2 c^TllV^^SI, (6.31) 

i=0 i=l 

where the constants Cj, c[ only depend on m and k. From this formula, we inductively obtain 
the desired bounds for \V k g\g. Similarly, the estimates for |V fc Rm(x,t)|g and \V k 4>(x : t)\g 
are obtained from (6.28) with a transformation analogous to (6.31). The lemma then follows 
by plugging everything together. □ 

Finally, we obtain our desired criterion for long-time existence. 
Theorem 6.13 

Let (<?(£), </>(i))te[o,T) solve (RH) a with non-increasing a(t) G [a, a], < a < a < oo and 
T < oo. Suppose that T < oo is maximally chosen, i.e. the solution cannot be extended 
beyond T in a smooth way. Then the curvature of (M,g(t)) has to become unbounded for 
t /*• T in the sense that 

limsup ( max|Rm(a;, t)\ 2 ) = oo. (6.32) 

t/-T VxGM / 
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Proof. The proof is by contradiction. Suppose that the curvature stays bounded on [0, T), 
say |Rm| < R . For any point x € M and vector X e T X M, define g(x, T)(X, X) := 
linit-^T g(x, t)(X, X). We estimate 

\g(x, T)(X, X) - g(x, t)(X, X)\ < J 2|S(x, r)(X, X)\dr < C\X\ 2 (T - t), 

where we used again the fact that § is uniformly bounded on M x [0,T). This shows that 
the limit g(x,T)(X, X) is well defined and continuous in x. Hence, we obtain a continuous 
limit g(-,T) € r(Sym 2 (T*M)) by polarization. From [4, Lemma 6.49], all metrics g(-,t) arc 
uniformly equivalent - as in (6.29) - which implies that this limit must be a (continuous) 
Riemannian metric. Moreover, we define (f>(x,T) := \im t ^T </>(£, 0j where we use again the 
embedding ejv : N R d to interpret as a map into M. d . We estimate 

\j>{x,T)-cj>{x,t)\< J \§- t <f>(x,T)\dT<C(T-t), 

since the bound on |V 2 </>| yields a bound on |Jj</>| = \t s 4>\. This implies that 4>{-,T) is well 
defined and continuous in x. The uniform bounds (6.27) from Lemma 6.12 then also hold 
for the limit (g(T),<p(T)) and hence g{T) and 4>(T) are smooth. Indeed, for an arbitrary 
background metric g, we have 

\V k g(T) - V k g(t)\ s < J \§- t V k g(r)\ s dr < C(T - t), 

which follows from the uniform bound for | JjVg|g that we have derived in the lemma above. 
This means, the convergence g(t) — > g(T) is smooth. With 

|V fc 0(T) - V k <f>(t)\s < [\it V k t(T)\ 5 dT < C(T - t) 

we see that also (f>(t) —> 4>{T) uniformly in any C fc -norm. Finally, restarting the flow with 
(g(T), 4>{T)) as new initial data, we obtain a solution (g(t), 0(t)) tS [r r +e ) by the short-time 
existence result from Chapter 2. This yields an extension of our solution beyond time T 
which is smooth in space for each time. From the flow equations and the uniform bounds 
on |V fc Rm| as well as |V fc ^>|, the time derivatives (and hence also the mixed derivatives) are 
smooth too, in particular near t — T . This means, the extension of the flow is smooth in 
space and time, contradicting the maximality of T. □ 

7 Monotonicity formula and no breathers theorem 

The entropy functional W Q introduced in this section is the analogue of Perelman's shrinkcr 
entropy for the Ricci flow from [31, Section 3]. It is obtained from the energy functional J a 
from (3.1), by introducing a positive scale factor r (later interpreted as a backwards time) 
and some correction terms. For detailed explanations of Perelman's result, we again refer to 
Chow et al. [5, Chapter 6] and Miiller [27, Chapter 3]. Moreover for the special case N C R, 
the entropy functional W a can be found in List's dissertation [22]. 
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The entropy functional and its first variation 

Let again g = g l3 6 T (Sym 2 + (T* M)) , <f> £ C°°(M,N), f : M -> R and r > 0. For 
time- independent coupling constant a(t) = a > 0, we set 



a07, 0, /, t) := y (r(i? s + | V/| 2 - a |V^) + / - m) (4iTr)- m ^e-f dV g . (7.1) 

As in Section 3, we take variations g e = g + eft., f £ = f + si, 4> £ = t^n{4> + such that 
<5g = ft, Sf = I and 8(f) = d. Additionally, set t £ = r + ea for some a £ R, i.e. <5r = cr. The 
variation 

SW a :=SW atg ,4j, T (h,4,i,o) := W a (<? + e/i,7rjv(0 + e0), / + ei,r + ea) 

CLE e—0 

is easiest computed using the variation of 3^ a and 



W a ( 5! ^/,r) = (H- m / 2 (rJ a ( 9 i,/)+ {f-m)e-Uv). (7.2) 

Using 

*/ (f-m)e-fdV=[ U+{f -m){\tr g h~ l)\e- f dV 

JM JM ^ ' 

and equation (3.2) for the variation of 3 a (g, <j>, f), we get from (7.2) 
5W a = / -rhij (Rij + ViVjf - aVi&Vjtydfi 

JM 

+ f r(| tr s h-t) (2A/ - |V/| 2 + i? - a|V0| 2 + ^)dfi 

JM 

+ [ (i + a (l- f )(R + |V/| 2 - a|V0| 2 ) - ^(/ - ro))<fc 

JM 

+ / 2ra 1 ?(r 9 ^-(V0,V/))d/i, 

where we used the abbreviation d/i := (4nT)~ m / 2 e~f dV . Rearranging the terms, writing 

£ = {-rh l3 + a 9ij )(^ gij ) + r(i tr 3 ft - £ - ^) (^), 
and using J M (Af - \V f\ 2 )d^i = -(Attt)-" 1 / 2 J m A( e -f)dV = 0, we get 

? a = / {-Thij + rrg, , ) ( /,', , + ViVj/ - CtVi^Vj^ - ^gij)dfJ, 
JM 

+ f t(\ tr fl ft - * - ^) (2A/ - |V/| 2 + i? - «|V0| 2 + t^)d[i 



5W 



M 



2Tad{ Tg 4>- (V0,V/))d/x. 



Fixing the background measure 

Similar to Section 3, we now fix the measure rf/i = (4 7 TT)- m / 2 e- / dU. This means, we have 
/ = - log((47rr) m / 2 !^) and from = 5dfx = (± tr g h-£-^)d/j,, we deduce I = \ tr s ft-ff . 
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Moreover, we require the variation of r to satisfy 5t = a = — 1. This allows us to interpret 
r as backwards time later. We then write 



■W£(g, </>, r) := W a (g, 0, - log {{^r) m ' 2 , r) 



(7.3) 



and 



<5K S At(^ *) : = * W a, ff ,*,- io g ((4. T )W^) ir (ft, | tr 9 ft + -1) . 
The variation formula above reduces to 

+ / 2TCM?(r s 0-(V0,V/>)<iM ! 
which is monotone under the gradient-like system of evolution equations given by 

f t g rj = -2(i2y + ViV 3 -/ - 
f0 = r 9 0-(V0,V/), 
^/ = -i?-A/ + a|V0| 2 + £, 



(7.4) 



(7.5) 



9< 
9< 



1. 





= -2Rc + 2aV</> <§ 


)V<f>, 




= ^ 




' — 7" 


= -1, 




, o 


= n*{(A-KT)- m ' 2 e 


- / ). 



As in Section 3, pulling back the solutions of (7.5) with the family of diffeomorphisms 
generated by V/, we get a solution of 



(7.6) 



Since W a is diffeomorphism invariant, we find the analogue to Proposition 3.1. 
Proposition 7.1 

Let (<?(£), 0(i))te[o,T) be a solution of (RH) a with a(t) = a > 0, r a backwards time with 

JLr 
dt ' 

the entropy functional W a (g, <p, /, t) is non- decreasing with 
d 

W a = 



-1 and solution of the adjoint heat equation under the flow. Then 



dt 



I 2t |Rc - a\7cj) ® V0 + Hcss(/) - (47rr)- m/2 e- / 
Jm 



+ / 2™!^^- (V0,V/)| 2 (47TT)- m / 2 e- / dV r . 



(7.7) 



Remark. As seen in Corollary 3.2 for the energy 1 a , the monotonicity of the entropy W a also 
holds true for non- increasing positive coupling functions a(t) instead of a constant a > 0. 

Minimizing over all probability measures 

Similar to \ a (g,(j>) defined in (3.9), we set 

fx a (g,c/>,T) :=inf {W£( ff ,0,r) | M (M) = 1} 



inf ^ W Q (g>, /,r) 



(47rr) _m/2 e _/ dV = 1 



,1/ 



(7.8) 
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Our goal is again to show that the infimum is always achieved. Note that for g = rg we 
have R S = ±R g , |V/|? = 1|V/||, |V0|? = \\Vcj>\% dV- g = T m / 2 dV g and thus 

lx a {rg,^,T) = n a (g,cj),l). 

We can hence reduce the problem to the special case where r = 1. Set v = (47r)~ m / 4 e~-^/ 2 . 
This yields 

W a (g,(f>,v,l) = J v(Rv - 4Av - a\V(j)\ 2 v - 2v\ogv - log(47r) - mvj dV, 

and hence fj. a (g, 0, 1) = inf{W Q (g, 0, i>, 1) | J M v 2 dV = 1} is the smallest eigenvalue of 
L(v) = — 4Av + (R— a\ V0| 2 — ^ log(47r) — m)v — 2v \ogv and v is a corresponding normalized 
eigenvector. As in Section 3, a unique smooth positive normalized eigenvector v m i n exists 
(cf. Rothaus [36] or List [22]) and we get the following. 

Proposition 7.2 

Let (g(t),<j)(t)) solve (RH) a for a constant a > and let J^r = — 1. Then fJ- a (g,(f>,T) is 
monotone non- decreasing in time. Moreover, it is constant if and only if 



= Rc - aV</> <g> + Hess(/) - -2 
O = r g </>-(V0,V/). 



2T ' (7.9) 



/or i/ie jnifiifTiizcT f that corresponds to Vmin . As always, the monotonicity result stays true 
if we allow a(t) to be a positive non-increasing function instead of a constant. 

Proof. The proof is completely analogous to the proof of Proposition 3.3, using the mono- 
tonicity of W„ instead on the monotonicity of 7 a . □ 

Non-existence of nontrivial breathers 

Breathers correspond to periodic solutions modulo diffcomorphisms and scaling, generalizing 
the notion of solitons defined in Definition 2.1. 

Definition 7.3 

A solution (g(t), </>(i))te[o,T) of (RH) a is called a breather if there exists t\,ti € [0, T), 
ti <.t<i, a diffeomorphism ip : M — > M and a constant c € R + such that 

f gfo) = cPg(t 1 ), 

\ 0(t 2 ) = 1 ' ' 

The cases c < I, c = 1 and c > 1 correspond to shrinking, steady and expanding breathers. 

Theorem 7.4 

Let M and N be closed and let (g(t), </>(i))te[o,T) be a solution of (RH) a with a(t) = a. 

i) If this solution is a steady breather, then it necessarily is a steady gradient soliton. 
Moreover, (j>(t) is harmonic and Rc = aV<p (8) V</>, i.e. the solution is stationary. 

ii) If the solution is an expanding breather, then it necessarily is an expanding gradient 
soliton. Again <p(t) must be harmonic (and thus stationary, <f>{t) = (f>(0)), while g(t) 
changes only by scaling. 
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in) If the solution is a shrinking breather, then it has to be a shrinking gradient soliton. 



If we assume in addition that dim M = 2 or that (M,g(0)) is Einstein, then in the first 
two cases above, <f>{t) is not only harmonic but also conformal, hence a minimal branched 
immersion, provided that it is non-constant. 



Proof. This is an application of the monotonicity results for X a {g,4>) from Proposition 3.3 
and for fi a (g,(j),T) from Proposition 7.2. Since the proof is very similar to the Ricci flow 
case solved by Perelman in [31], we closely follow the notes from Kleiner and Lott [19] on 
Perelman's paper. 

i) Assume (g(t), ^>(i))te[o,T) is a steady breather. Then there exist two times t\, ti 
such that (7.10) holds with c = 1. From diffcomorphism invariance of X a (g,(f) 
defined in (3.9), we obtain X a {g, <fi)(ti) — X a (g,<p)(t2). From Proposition 3.3, we 
get condition (3.10) on [ii,^], which means that (g(t), <fi(t)) must be a gradient 
steady soliton according to Lemma 2.2 and uniqueness of solutions. Moreover, the 
minimizcr / = — 21ogu rra i„ which realizes X a (g,4>) is the soliton potential. From 
(-4A + R - a\V 4>\ 2 )v rnin = X a (g, 4>)v min =: X a v min , we obtain 

2A/-|V/| 2 + i?-a|V0| 2 = X a . (7.11) 

Since (g(t),4>(t)) is a steady soliton, (2.5) holds with a = 0. Plugging this into (7.11) 
yields A/ — |V/| 2 = X a , and we obtain from J M e~^dV = 1 

X a = f X a e~ f dV= f (A/-|V./f)e-^y = - f A(e~ f )dV = 0, 

JM JM JM 

i.e. A/ = |V/| 2 . Another integration yields 



/ \Wf\ 2 dV= f AfdV = 0, 
Jm Jm 



' M JM 

and thus V/ = 0, Hess(J) = on M x [0,T) and (3.10) becomes 

J = Rc - aV(j)® V0, 
\ = r g 4>. 

In particular, </>(i) is harmonic and (g(t), 4>(t)) is stationary. 

ii) The proof here is analogous to the case of steady breathers, but we first need to 
construct a scaling invariant version of X a (g, 0). We define 

/ r \2/m 

X a {g,4>):=X a {g^)[JdV g \ . (7.12) 

This quantity is invariant under rescaling g = eg. A proof of this fact is given in the 
appendix of the author's dissertation [28]. Moreover, we claim that at times where 
X a (t) := X a (g,(j))(t) < 0, we have JjA Q (i) > 0. Indeed, note that / = \ogV(t) satisfies 
J M e~fdV = 1 and is thus an admissible test function in the definition of A Q , hence 

\ a (g,<t>)<? a (g,cj>,logV(t))= [ Se-^^dV^Vitr 1 f S dV. (7.13) 

Jm Jm 
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With the assumption \ a (t) < 0, we find 

& X* > 2V 2 ' m j (|S + Hess(/) - i(5 + A/) ff | 2 + «|r g( /» - (V<£, V/)| 2 )d/z 



(7.14) 



the right hand side being nonnegative by Holder's inequality (see again [28] for a more 
detailed computation). Now, assume that (g(t),<j){t)) is an expanding breather. Since 
A a (g,4>) is invariant under diffcomorphisms and scaling, we have A a (tx) = Xafo) for 
the two times t\, t% that satisfy (7.10). Since V(h) < V(t 2 ), there must be a time 



to e [h,t 2 ] with -§jV(t ) > and hence with (7.13) 



Aa(*o) < Vito)- 1 / S dV = -Vito)- 1 ^^) < 0. 

JM 

The claim applies and we obtain X a (ti) < X a (to) < and since X a (t 2 ^ — A Q (ti), we 
see that X a (t) must be a negative constant. Hence, both lines on the right hand side 
of (7.14) have to vanish. This means that (S + A/) has to be constant in space for all 
t and because X a (t) = J M (S + A/)d/x, this constant has to be X a (t). From the first 
line of (7.14) we obtain 

fo = Rc-«V</)®V0 + Hess(/)-^ 3 , 
\O = T^-(V0,V/). 

By Lemma 2.2, (g(t), <^(i))te[o,T) is an expanding soliton with potential / = — 21ogu m m- 
This means that we can use (7.11), which implies 

= 2A/ - |V/| +S-X a = 2A/ - |V/| + S - (A/ + S) = Af - |V/| (7.16) 

and thus by integration V/ = 0, Hess(/) = 0, as above. Plugging this into (7.15), the 
second equation tells us that <fi{t) is harmonic and the first equation yields 

§- t g = -2Rc + 2aV0 ® V0 = -2%g, 

i.e. (M,g(t)) simply expands without changing its shape. 

iii) If (g(t), 4>(t)) is a shrinking breather, there exist t\, t 2 and c < 1 which such that (7.10) 
is satisfied. We define 

t 2 — ct\ , . 

t := — > t 2 , and r(t) = r — t. 

1 — c 

Note that r(t) is always positive on [£i,i2]- Moreover, c = (to — t 2 )/(ro — <i) = 
t(^2)/t(£i). Then, from the scaling behavior of /i Q (<7, 4>, t) and diffcomorphism invari- 
ance we obtain 

Va{g(h), 0(*2), t(* 2 )) = fj, a (cip*g(tx), tp*^(h), ct(*i)) 

= H a (ip*9(h),ip*<l>(ti),T(ti)) (7.17) 
= /^(g(£i),<?Kii),T(ii)). 

By the equality case of the monotonicity result in Proposition 7.2, (g(t),4>(t)) must 
satisfy (7.9) and according to Lemma 2.2 thus has to be a gradient shrinking soliton. 
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It remains to prove the additional statement in the cases where dimM = 2 or (M,g(0)) is 
Einstein. If (g(t),(f>(t)) is a steady or expanding breather, we have seen that -^g — eg. In 
particular, if (M,g(t)) is Einstein at t = it remains Einstein under the flow. Moreover, 
since Rc = — q in these two cases, we get 

Wl)n = V^V^ = + 2%) = ^(2£ + c ) 5y , 

i.e. is conformal. It is a well-known fact that conformal harmonic maps have to be minimal 
branched immersions (cf. Hartman-Wintner [16]). □ 

8 Reduced volume and non-collapsing theorem 

Let us briefly restate the main result from our previous article [29] about the monotonicity 
of reduced volumes for flows of the form j^gij = —2Sij, where Sij is a symmetric tensor 
with trace S = g^ Sij. (The resuts can also be found in the authors thesis [28].) 

Monotonicity of backwards reduced volume 

In order to define the backwards reduced distance and volume, we need a backwards time 
r{t) with Jjr(i) = — 1. Without loss of generality, one may assume (possibly after a time 
shift) that t = —t. 

Definition 8.1 

Assume -jfpgij = 2Sij has a solution for r G [0, f] and < t\ < < f, we define the 
Lb-length of a curve r\ : [tq,ti] — > M by 

MV) ■= r (S(v(r)) + \£v(t)\ 2 ) dr. 

J To 

Fix the point p G M and Tq = and define the backwards reduced distance by 



£b{q,Ti) := inf 
776 r 



where T = {77 : [0,ri] — > M \ 77(0) — p, 77(71) = q}. The backwards reduced volume is defined 
by 

Vf,(r) := f {AKT)- m/2 e~ lbi - q ' T) dV{q). (8.2) 

J M 

The following is proved in [29, Theorem 1.4]. 
Theorem 8.2 

Suppose that g(t) evolves by ^g%j = —2Sij and the quantity 



V{S,X) := -j^S — AS — 2 |Syf + 4(ViS tf )*i - 2(V j S)X j 
+ 2RijX{Xj — 2SijX{Xj , 



(8.3) 



is nonnegative for all vector fields X G T(TM) and all times t for which the flow exists. 
Then the backwards reduced volume V^(r) is non-increasing in t, i.e. non- decreasing in t. 
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Ill our case where Sij is given by Rij — ctVj^V j<j>, the evolution equation (4.14) for Sij 
together with ^{V l S lj )X j - 2(VjS)X,- = -AaT g 4>Vj<j>Xj yields 

£>(,%, X) = 2o\T g 4> - V X (j)\ 2 - a\Vcj)\ 2 

for all X on M. This means, T>(Sij,X) > is satisfied for the (RH) a flow with a positive 
non-increasing coupling function a(t) and the monotonicity of the reduced volume holds. 

No local collapsing theorem 

We have seen in Section 6 that the metrics g(t) along the (RH) a flow are uniformly equivalent 
as long as the curvature on M stays uniformly bounded. But a-priori, it could happen that 
at a singularity (i.e. when Rm blows up) the solution collapses geometrically in the following 
sense. 

Definition 8.3 

Let (g(t), (f>(t))te[0,T) °e a maximal solution of (RH) a , or more generally of any flow of the 
form -j^gij = —2Sij. We say that this solution is locally collapsing at time T, if there is a 
sequence of times tk /* T and a sequence of balls Bk ■= B g u k \(xk,rk) at time tk, such that 
the following holds. The ratio rt/tk is bounded, the curvature satisfies |Rm| < r^ 2 on the 
parabolic neighborhood Bk X [tk — and r^ m vol(Bk) — > as k — > oo. 

Using the monotonicity of the reduced volume, we obtain the following result. 
Theorem 8.4 

Let (g(t),(j>(t)) be a solution of {RH) a with non-increasing a(t) £ [a, a], < a < a < oo 
on a finite time interval [0, T). Then this solution is not locally collapsing at T. 

The only ingredients of the proof are the interior gradient estimates from Corollary 6.11 
and the monotonicity of the backwards reduced volume stated above. Hence, every flow 
Jj<7 = — 2§ that satisfies the assumption of Theorem 8.2 and some interior estimates for S, 
VS in the spirit of Corollary 6.11 will also satisfy the non-collapsing result. For the (RH) a 
flow, it is possible to obtain a slightly stronger result using the monotonicity of /i Q (g, 0, r) 
from Section 7 instead of the monotonicity of the backwards reduced volume. In the special 
case N C R, this can be found in List's dissertation [22, Section 7]. The proof in the case 
of (RH) a is analogous. However, the result here is more general in the sense that it may be 
adopted to other flows -^g = — 2§ in the way explained above. 

Proof. The proof follows Pcrelman's results for the Ricci flow [31] very closely, see also the 
notes on his paper by Kleiner and Lott [19] and the book by Morgan and Tian [25]. However, 
we need the more general results from [29] that also hold for our coupled flow system. We 
only give a sketch. 

The proof is by contradiction. Assume that there is some sequence of times tk T and 
some sequence of balls Bk := B g ^ tk ^{xk,rk) at each time tk, such that r\ is bounded, the 
curvature is bounded by |Rm| < r^ 2 on the parabolic neighborhood Bk x [tk — r \^k\ an d 
r A 7 m vol(i?fe) — > as k — > oo. Define Sk ■= r^ 1 vol^fc) 1 /™, then e fc -> for k — > oo. For 
each k, we set Tk(t) = tk — t and define the backwards reduced volume Vk using curves going 
backward in real time from the base point (xk, tk), i.e. forward in time Tk from Tk = 0. The 
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goal is to estimate the reduced volumes Vfc(efe^), where rfc = Skr\ corresponds to the real 
time t = tk — £tr'l, which is very close to tk and hence close to T. 



Claim 1: limfe-^ Vk{e k r 2 k ) = 0. 



Proof. An £;,-gcodesic tj(t) starting at 77(0) = Xk is uniquely defined through its initial 

8 fc fe 



vector v = lim T _j.o 2^/tX = \m\\^ a X. First, we show that if \v\ < \£ k X ^ 2 with respect 



1 II 

to the metric at (xk,tk), then t](t) does not escape from B k := B g u k )( x k,'i"k/2) in time 
t = e k r\. Write ik — t k — r\. Since |Rm| < r k 2 on B k x by assumption, we obtain 

from Corollary 6.11 

|V0| 2 < — |V 2 0| 2 <-^-, |VRm| 2 <— ^— , on Bl /2 x(i k ,tk), 

t-tk [t-t k y {t-t k y 

for some constant C independent of k. Without loss of generality, Sk < \ so that t — tk > \r\ 
whenever t £ [t k — E k rf,,t k ). This means that 

|V0| 2 <Cr k 2 , \V 2 (b\<Cr k 2 , |VRm| < Cr k \ on b\' 2 x [t k - e k r 2 k , t k ). 

Together with the assumption |Rm| <r k 2 , this yields 

[SI < |Rc| + IV0I 2 < CrT 2 , 
|V^| < |V,R| + |V0||VV| < C 1 ^ 3 , 

on -B^ 2 x [tk — £ k rl,t k ). Plugging this into the estimate (4.5) of [29], we get 

JL\X\ < \C\X\r k 2 + X 2 Cr k 3 < C\X\e]j 2 r k l + Ce k r k \ (8.5) 

for A = y/r < \fekTk = £ k 2 ' rk - Since |X(0)| = \v\ < we obtain the estimate 

\X{\)\ < i£fc 1/2 for all r G [0,e fc r 2 ] if k is large enough, i.e. Sk small enough. With an 
integration, we find 

/ \X(r)\dr = \X(X)\d\< l e -^ 2 d x<\ rk . 

Jo Jo Jo 

Since the metrics g(r = 0) and g(r = £kr 2 ) are close to each other, the length of the curve 
77 measured with respect to g(r = 0) = g(tk) will be at most r k /2 for large enough k. This 
means that indeed 

(j?(T),t fc -t) G B^ /2 x [t fc - efcrgjtfc), V0 < t < e fe r^. (8.6) 
With the bounds from (8.4) and the lower bound in [29, Lemma 4.1], we obtain 
Mv) > ~Cr k 2 (e k r 2 k )V 2 = -Cel /2 r k , i.e. £ b (g, e k r 2 k ) > -Cs k 



Thus, the contribution to the reduced volume ^(e^r 2 ) coming from £b-geodesics with initial 
vector |u| < ^e k ^ 2 is bounded above for large k by 



/ {^e k rl)- m/2 e Cek dV < Ce k m/2 r k m vol(B, 1/2 ) < Ce k /2 -> (k -> 00). 
Jb] /2 
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Next, we estimate the contribution of geodesies with large initial vector \v\ > g£fc to 
the reduced volume V k (e k r k ). Note that we can write the reduced volume with base point 
(xk,tk) as 

V k (n) = [ (47TTi)- m/2 e-^' T1 W(<7) = / (47rri)- ro/2 e-^ 6exp ^ W ' Tl) J(i;,T 1 )d'i;. 
Jm Jn{T!,k) 

Here, X^exp^ is the ^-exponential map defined in [29], taking v to 77(7-1) with 77 being the 
£b-geodesic with initial vector v, J(v, n) = det d(£&exp" ) denotes the Jacobian of -Cfeexppj, 
and f2(ri, fc) C T Xk M is a set which is mapped bijectivcly to M up to a set of measure zero 
under the map i^exp^. In [28], we prove that the integrand 

f(v,n) := (47rri)- m / 2 e - £(£60xp ^ (u) ' ri) J(7;,Ti) (8.7) 

is non-increasing in n for fixed v and has the limit lim Tl _>.o f(v, 7i) = 7r _, ™/ 2 e _ ' t '' . Together 
with f2(r') C f2( r ) for r < r' this yields an alternative proof of the monotonicity of the 
reduced volumes obtained in [29]. Moreover, it implies that the contribution to the reduced 
volume V k {£ k rj,) coming from £{,-geodesics with initial vector \v\ > |e fc can be bounded 

n -m/2 e -\v\' dv < Ce - 647I ^ (fc^oo), (8.8) 

M>K 1/2 

which completes the proof of Claim 1 . □ 
Claim 2: V k (t k ) is bounded below away from zero. 

Proof. Let us remark that r = tk corresponds to real time t = 0. We assume that k is large 
enough, so that tk > T/2. The idea behind the proof is to go from (xk,tk) to some point qk 
at the real time T/2 (i.e. r = — T/2) for which the reduced ^-distance £b(qk,tk — T/2) 
is small. From the upper bound on Lb from [29, Lemma 4.1], we see that for small r it is 
possible to find a point <7fc(r) such that ^(,(%(r), r) < ^. On the other hand, combining the 
evolution equations for Jp4 and A4, we obtain 

£| T=ri 4 + A4<-^4 + ^ (8.9) 

(in the barrier sense) and hence for the minimum of l m in(j) = niin ge M 4(<7, t) 

"Q^\ T — Tl ^min ~^min T gri (8.10) 

in the sense of difference quotients. The latter is obtained by applying the maximum prin- 
ciple to a smooth barrier. The inequality (8.10) shows that there is some point qkij) with 
£b(qk(T~), t) < y f° r every r. As mentioned above, we choose qk at the real time T/2 with 
4(?fcTfc - T/2) < y- Let 77 : [0,t k - T/2] -> M be an £ 6 -geodesic realizing this length. 
Moreover, let 77^ : [tk — T/2,tk] —> M be g(t = 0)-gcodcsics (i.e. a g(r = ife)-geodesic) from 
q k at time r = t k - T/2 tope B qk := B g{T=tk) (q k ,l) = B g{t=Q) (q k ,l) at time r = t k . 
Since |Rm| is uniformly bounded for t G [0, T/2], we get a uniform bound for S along this 
family of curves. Since all the metrics g(r) with r G [tk — T/2,t k ] are uniformly equivalent, 
we get an uniform upper bound for the £b-lcngth of all 77^. From this, we see that the 
concatenations (77 - — ' 77^,) : [0,tk] — > M connecting Xk to p G B qk have uniformly bounded 
£&-length, independent of p and k. This gives a uniform bound h(p, tk) < C, for all p G B qk 
and k G N large enough. We can then estimate 

V k (t k )= I (4TTt k y m / 2 e- lb{q < tk) dV(q)> [ (4TTt k y m / 2 e~ c dV >C inf vo\(B qk ), 
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which is bounded below away from zero, independently of k. This proves Claim 2. □ 

Since the backwards reduced volumes Vk are non-increasing in t (i.e. non-decreasing in real 
time t) according to Theorem 8.2, we obtain Vfc(ifc) < ^fe( £ fc r fc) f° r & large enough. But 
since Vfc(ifc) is bounded below away from zero by Claim 2 while Vk(£ktf.) converges to zero 
with k — > oo by Claim 1, we obtain the desired contradiction that proves the theorem. □ 

A Commutator identities 

It is well-known that a (p,q)-tensor B (i.e. a smooth section of the bundle (T*M)® P <£> 
(TM)® q ) satisfies the following commutator identity in local coordinates (x 1 ,...,x n ) in- 
duced by a chart : U -> M", U C M, 

q p 
[\7 V7 1 R fc i-" fc 9 _ r>k r u k 1 ...k r - 1 mk r+1 ...k q ST^ „ „ki...kg 

[Vi,VjJ£f £l ^ - 2_, H ijm ti ll l ...Z p +2 J %<.m%..i,_ 1 m<, +1 ..i p ' (A.i) 
r— 1 s— 1 

Now, assume that we are given Levi-Civita connections for all (p, g) tensors over (M, g) and 
over (N, 7). For a map : (M,g) — > (iV, 7), there is a canonical notion of pull-back bundle 
<f>*TN over M with sections 0*^ = Uo for V G T(TN). The Levi-Civita connection V TJV 
on TJV also induces a connection V* T7V on this pull-back bundle via 

v£ TJ VV = <f(V£^F), X G r(TM), V G r(TJV). 

We obtain connections on all product bundles over M with factors TM, T*M, cjfTN and 
(f)*T*N via the product rule and compatibility with contractions. Take coordinates x k on 
M, k — 1, ... ,771 = dimM, and y M on iV, /i = 1, . . . , n = dim N, and write dk for and 
d„ for jjL We get ViVjV K - V J V l V K = R tJKX V x with 

Ri jK x(x) = (Rm(d i ,d j )<l>*(d x ),cf>*(d K ))^ TN (x) 
= ( w Rm( ^*5i,</)* dj)d\,d K ) TN (0(x)) 
= \,kA V< ^ (s) Vj 0" (x) , 

where we used 0*9, = V^0^ <9 M . This allows to extend (A.l) to mixed tensors, for example 

[Vi, Vj]-B^ = Ri jp B^ + RijipB^ + Ri je B^ + RijxgB^. (A.2) 

The standard example that will be used quite often is the following. The derivative \7 <fi of 
: M — > TV is a section of T* M (g> 4>*TN . Thus, the intrinsic second order derivative is built 
with the connection on this bundle, i.e. V.;V,-c/> A = didj4> x - T^d k (j) X + ^^d^dj^ and 
similar for higher derivatives. Using (A.2), we obtain 

ViVjVi^ - V,V,;V £ / = RijgpVp^ + X.aV^V^V^. (A.3) 

There is also a different way to obtain these formulas, which is especially useful when </> is 
evolving and we also want to include time derivatives. We learned this from [20]. Here, 
we interpret \7 k <fi as a k- linear TN- valued map along <f> £ C°°(M,N) rather than as a 
section in (T*M)® fe ® <j>*TN. Leting ui be any such fc-linear TTV-valucd map along 0, 
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i.e. u>(x) : (T x M) xk — > T^^N, the covariant derivative Vw is a (k + l)-linear TN- valued 
map along (f> etc. The curvature tensor fc Rm for u> can then be computed by 

( fe Rm(A, Y)u)(X u ...,X k )= N Rm{Vct?{X), V#y))w(Xi, . ..,X k ) 

* (A.4) 
-^2u(X 1 ,...,-Rm{X,Y)X„...,X k ). 

s=l 

Of course, this agrees with the definition above, where we used the bundle interpretation. 
Now, if <j) is time-dependent, we simply interpret it as a map <fi : M x I — > N and interpret 
W k (f> as fc-linear TiV-valued maps on M x I along <f>. The formalism stays exactly the same. 

Note that induces a covariant time derivative Vt (on all bundles over M x I) that agrees 
with for time-dependent functions. Choose coordinates x l for M with 

V t (&) = Vi^-) = Vt(fls) = Vi(&) = 0, Vi,j = l,...,m (A.5) 
at some base point (p, i) in M x J. Then, using (A.4), we obtain for u> = 

Vf(ViV^) = Vt((ViCj)(^)) = (V t V^)(9,) = (V 1 V( W + 1 Rm(|,a i ) W )(9,) 

^,V^)o;(^)-a;( MxJ R m (^ 



Remark. If we also vary the metric g on M in time, we will get an additional term from the 
evolution of ViVj, namely — (^r^)Vfe0. Note that is a tensor, while T itself is not. 
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